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ABSTRACT , 

This thesis presents a procedural Pe cice for processing and 
assessing the integrity of experimental strain gage data taken on piping 
in which gage rosettes are arranged in circumferential patterns at 
several sections of a piping system. The method of least squares is 
used to obtain "best" values for strain per unit loading for each gage 
element. This data is used to find principal stresses and values of 
normal and shearing stress on a normal cross section at each rosette 
location. The data is then smoothed to obtain normal and shearing stress as 
functions of angular position for each ring of rosettes, A computer 
program has been devised to perform these operations and provide a 
graphical output to assess degree of fit with experimental data, Finally, 
manual computation is employed to study the equilibrium of various seg~ 
ments of the piping in evaluating the overall integrity of the data, 

This procedure and program have been used in evaluating tests performed 
on the main seawater piping system of the USS Benjamin Franklin (SSBN 
640). Some results of the dockside hydrostatic tests are presented 


herein for purposes of illustration. 
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1. Introduction. 
“ee .. ORIGIN OF IDEA FOR THESIS 

The idea for this thesis was conceived while the writer worked at the 
Klectric Boat Division of General Dynamics Corporation during an indus-. 
trial tour from the U. S. Naval Postgraduate School, 

While working on the analysis of strain gage data taken during docke 
side hydrostatic tests and sea trials of the USS Casimir Pulaski (SSBN 
633), it became evident that a method to fully analyze the main seawater 
piping system from the strain gage data did not exist. The only tool 
available at the time was a computer program that had been developed by 
Electric Boat which determined the two principal stresses and principal 
stress axes for each individual set of rosette readings. 

The project of instrumenting the main seawater piping system on the 
USS Benjamin Franklin (SSBN 640) and taking strain gage data during dock- 
side hydrostatic tests and sea trials had been assigned to Electric Boat 
by the U.S. Navy Bureau of Ships, but no assignment had been made to 
analyze the data taken, The writer obtained the necessary permission 
from Electric Boat and the Office of the U. 5. Navy Supervisor of Ship- 
building, Groton, Connecticut, to use the data in connection with a thesis 
at the U. 5. Naval Postgraduate School and established liaison with 
personne] at Electric Boat connected with conducting the SSBN 640 tests. 

Over the last nine months Mr, Carl P. Wilson, Supervisor of the 
piping section of the Marine Engineering Design Department at Electric 
Boat and his main assistant, Mr. William R, Hintz, have been most helpful 
and prompt in forwarding the test data and other related material requested, 
Without their assistance and cooperation this thesis could not have been 
possible. The writer would like to express his sincere appreciation to 


Dr. John &, Brock, Professor of Mechanical Engineering, U. S. Naval 
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Postgraduate School, for his continued patience, encouragement, and most 
capable guidance while acting as faculty advisor, 

The intention of this thesis is not only to analyze the data from the 
SSBN 640, but also to develop a procedural sequence for processing and 
assessing the integrity of strain gage data that could be utilized in the 
design and testing of future submarine seawater piping systems. 

. HOW DATA WAS OBTAINED 

Two complete sets of data were used jim developing this thesis. 

The first set of data was taken on 7 Aygust, 1965 during the dockside 
hydrostatic test of the main seawater piping system of SSBN 640 in 
accordance with Electric Boat Division's Propulsion Plant Test Form 

Number 721-5 Revision A, The second set of data was taken during builder's 
sea trials on 19 and 20 September, 1965 in accordance with Electric Boat 
Division's Sea Trial Form Number 21.04 Revision A. However, for security 
reasons, no further reference is made in this theses to sea trial data. 

All strain gages utilized were BLH Type FABR-25-1259, constantan grid, 
phenolic glass base, three element 45° superimposed grid rosettes. More 
detailed information concerning these rosettes may be obtained from BLH 
Specifications Catalog Bulletin 101 effective September 1964. 

Gages were equally spaced around a cross section of piping at 30° 
from each other, A cross section with twelve rosettes will be referred to 
as a ring or rosette ring throughout the remainder of this thesis. The 
rosettes in each ring were oriented as shown in Figure 1.2.le,, There were 


eight such rings. 
FLOW 
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The number twelve gage is located facing directly up for all rings 
except ring number <, in which it faces port. The approximate locations 
of the rosette rings are shown in figures 1.2.2, 1.2.3 and 1.2.4 which are 
sketches of the piping system. The exact location of each ring and the 
piping system dimensions are not given, The dimensions may be found on 
Bureau of Ships Blueprints SSBN 640 209 2117801 Rev. H and SSBN 640 209 
2117802 Rev. F and the exact ring locations may be obtained from either 
Electric Boat Division's Propulsion Plant Test Form Number 721-5 Rev. A 
or Sea Trial Data Form Number 21.04 Rev. A. 

Strain gage readings were taken using SR-4, Strain Indicators Type N 
manufactured by Baldwin-Lima-Hamilton Corporation. 

OBJECTIVES OF THESIS 

The basic purpose of this thesis is to develop methods of assessing 
the integrity of the data and, incidentally, to reduce the data to a form 
more convenient for subsequent analysis, To reach this goal a digital 
computer program has been developed and papleged. This program uses as 
imputs certain constants concerning the piping material and sizes and 
the raw strain gage data, It computes the two principal stresses, the 
angle measured from the direction of flow through the piping of the 
maximun principal stress, and the shear stress at each cross section and 
stress normal to each cross section. In addition, smoothed values of 
normal and shear stress every..OQl radian around the circumference of the 
rosette rings are calculated and plots are drawn showing both smoothed and 
actual data points, Finally the program computes the three forces and 
moments at each rosette ring. 

The outputs of the computer program are analyzed to determine the 
integrity of the raw data, if the maximim principal stress exceeds yield 


strength, the agreement between smoothed and actual data values of 
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normal and shear stress and if equilibrium between two rosette rings is 
established, 

For the following sections, the terms experimental data and 
experimental results are used to represent the actual strain gage data 
taken and the computed results from this data, 

Actual dimensions of the piping system, data, and actual numerical 
results have been eliminated from this thesis in order to keep the 
thesis unclassified. All mumerical results given are the actual results 
multiplied by arbitrary constants. However, relative values of the 
results are the same as the actual results. 

: NOMENCLATURE 


a,b,c,d, } Constants for smoothed value of normal stress 
e,f,g 


E Young's Modulus 
ei Error 
Cr , Ce Unit vectors in the 7 and ~ directions 


F(¢:),F,(@:) Functions of the angle g 


Fax 

Fry Resultant equilibrium forces in the x,y and z directions 

F Re 

Fy Force on a cross section of the main seawater piping system 
in the direction of the ship's bow, 

Fy Force on a cross section of the main seawater piping system 
in the vertical direction from ship's centerline. 

Fa Force on a cross section of the main seawater piping system 
in the starboard direction from ship's centerline. 

Fy 

Fy! { Forces in the x’ , y’ , and z’ directions 

F 3! 

i,j,k Unit vectors in the x, y and z directions. 


k,l,m,n, } Constants for smoothed value of shearing stress 
OsP39 


X,y¥,23X',y4 z/ Coordinate systems described in text 
i 


Mrx 
Mery 
Mere 


Mx 


y 
Lt J 
[a] 


(B] , (B'] 
(ce), (Cc 
Son (Se E3 
Op 

p 


ys 


Imax 


Omin 


On 


iM? ¢ 


Resultant equilibrium moments in the x, y and z directions, 
Moment about the fore and aft axis of the ship in the 
direction of the bow, 


Moment about the vertical axis of the ship in the upward 
direction. 


Moment about the port and starboard axis of the ship in the 
direction of starboard 


Moments in the x' , y’ , and z! directions. 


Radius of the piping. 

Pipe wall thickness, 
percent of test depth. 
Strain reading at any depth 
Matrix 

a 7x 7 matrix. 


7 x 1 Matrices 


Three strains of a three element 90° rosette. 
Maximum principal stress angle 

Poisson's ratio. 

Summation sign 

Maximum principal stress 

Minimum principal stress 


Stress normal to a cross section of piping in the direction 
of filows 


Shear stress on the outside surface of the pipe positive in 
the direction of increasing rosette numbers, 


Shear stress in the xz plane 
Angle on piping cross section from rosette number 12, 


Stress tensor (dyadic). - 


n2 


2, Theory! 
THEORY FOR STRAIGHT LINE FITTING 
It is assuned that the strain gage readings will increase linearly 
with internal pressure which in turn does increase linearly with an in- 
crease in depth. This assumption is based on the following facts: 

(a) The piping material exhibits a straight line stress-strain 
Curve. 

(b). The system is not loaded above the elastic limit of the 
material as verified by the maximum values of the maximum principal stress 
at every rosette. 

(c) The deflections are small as verified by the overall change 
in strain gage readings. 

(d) It is assumed, for the builder's sea trials, that hull 
deflection inputs to the piping increase linearly with depth. For dock- 
side hydrostatic tests, there was no such input. 

The other main effect acting on the system is the weight of the system. 
Since the system is always filled with seawater of nearly constant density 
the weight effect on the piping system will remain nearly constant and, 
therefore, not have any appreciable effect on the assumption of linearity. 

The method of least squares is used to obtain the best straight line 

connecting all data points of each element with the percentage of test 
depth taken as the independent variable and the strain gage readings 

taken as the dependent variable. The results of the method of least 


squares are given in the form of the equation: 


Ioreater details of the theoretical relations are given in, or may 
easily be inferred from, corresponding subsection of section 3, which deals 
with the computer program used to perform the calculations, 


: 2spiegel, M.R., Schaum's Outline Series, Theory and Problems of 
Statistics (New York; Schaum.Publishing Co., 1961) p.220 
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Y a + bx 2c la 
where a is the zero percent test depth intercept and b is the slope of 
the best straight line connecting all data points from the surface to 
test depth. The slope, b, is also proportional to change in strain per 
percent changes in test depth. The slopes are the values of the indi- 
vidual element strains, @), G2, and &3, used in the subsequent develop- 
ment of the theory. 

THEORY FOR SOLVING PRINCIPAL STRESSES AND ANGLES 
The equations defining the principal stresses in terms of the 
strains measured by 4 rectangular rosette ares 


max = Ef Enea + We, €.)- Hote “Ce, ee neal 
2 


1—- AL a 


Omin = §} £6 — ley (€,-€2)" + [2€.-(e.+ €s)]* | 2. ome 


dmAL Leh 


where the elements of the rosette are numbered as shown in figure 2.2.1. 





2 





as — Direction of algebraically 
greater stress. 
©p 


Figure 2.2.1 
The direction of the maximum principal stress from element number one 


of the rosette in terms of the measured strains is: 


Gp= Larctan 2€2 -(€,+6€3) roe 
e €,- €3 


The derivation of equations 2,2.1, 2.2.2, and 2.2.3 is given 
in Perry and Lissner, The Strain Gage Primer, pages 136-138 


Li, 


’ THEORY FOR EXPERIMENTAL VALUES OF Ow # Txy 
Looking down on the surface of the piping, the stress picture looks 
like figure 2.3.1. 


direction of flow 





Omin 


3e— «w- 


Om Ax 


Figure 2.3.1 
It is known that G3 which is perpendicular to the plane of O and G2 
is equal to zero. 
To determine the stresses on a plane perpendicular to the x-y plane 
whose trace is mn, tensor theory is used, Temporarily use coordinates 
€, 1, and z to determine O- N., tie and tx2. (The ence 2 is per- 


pendicular to the plane shown in figure 2.3.1 with unit vector k.) 


From tensor theory, the stress tensor (dyadic) may be written: 


ef = Tinax Egle + Grn Cy ey ~ oem 
Gre beget BSD 


Tay = beds F VAS & 
tusk ge 


Since the term Cif appears in equations 2.3.2, 2.3.3, and 2.3.4, solve 
for it first. 
Ted = Cmax (C- Ce) €e + Gmin (¢-€7) ex 

= Umax Cos Op es — Umin Sin Op Cr Viggo. 
Substituting equation 2.3.5 back into equations 2.3.2, 2.3.3, and 2.3.4 
results in the following equations: 

Ou= Omax Cos Gp (Ce i) - Gmin SIN Op (7: i) 
= Omnax cos” Op + Gmin Sin*-Op Zoe 
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oe = Omax COS Op( Ce -j) - Omin SIN @p (Ex-j) 
= Omax COSOp SING p - Umin SiIN@p COS Op 


/2 (Gmax - Tm) SIN 28p 2.3.0 
Tx2z = Tmax Cos@p (€e-k) -~Gmin SInGp (€n-k) 
(Gmax COS Op) O -(Gmw SIN@p) O 


O 2B 


be ~-PREORY FOR SMOOTHED VAIUBS OF NORMAL AND-GHEARING STRESS 

In developing the theory for the smoothed values of Gnand Txy 
it has been assumed that the radius of the piping has remained 
constant throughout the entire piping system. This assumption is based 
on the diameter readings taken by GD/EB of rosette rings two through 
eight. These readings indicated that the diameter of the piping around 
a rosette ring remained constant enough to use a constant value of radius 
for calculations in this thesis. 

Since a previously developed theory was not available, it seemed 


reasonable to assume that Oy and Txy vary only with g in the following 


Fourier Series; 


In= A+ bsinG+ ccos@ +A SiIN2D +ecos 29+ sin 3D 


+GCos BG@+ «+++ = F(Q@) 2 ae 
Txy= k+ Ising@+ mMcosfZ+NSINZQ+ OC COSzP+ PSNI 
+ GCOS BG +111) = F,\(D) 2.4.2 


One of the decisions that had to be made concerned the number of 
terms of the Fourier Series to use. Since curves showing the exact 
shape of the stress distribution around the circumference did not exist, 
there was no criterion available to which curves developed could be 
compared, To assist in deciding how many terns to use, the solution was 


investigated using five, seven, nine, and eleven terms, A plot of the 
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experimental values and the smoothed values of Jy for five, seven, nine, 
and eleven terms for rosette ring number three is shown in figure 2.4.1. 
A decision was made to use seven terms for the following reasons: 

(a) A certain — of redundancy was necessary to permit 
good least squares evaluation. 

(b)If higher soperi clenesavare really of significance, more 
than twelve rosettes would be required at each rosette ring to obtain 
more than seven coefficients whose accuracy could not be doubted. 

(c) For equilibrium, only the first three coefficients are used, 
Therefore, it is necessary to have these evaluated as well as possible. 
The truncated forms of equations 2.4.1 and 2.4.2 shown below will be 
used in the subsequent development of theory. 


One 2+ bsngd+ceosft+Asin2hG+ecos2d + Fsin3¢d 


+9 cos 3G = F(@) 2,4 ae 
Try = k+ Lsing +mMeoshG+NSIn2P + OCOSZG+ PSIN 3 
+qGcos3g = Fi(Q) 2.4.2@ 


The smoothed values of Gx will be discussed before proceeding with the 
smoothed values of Txy. 

There are N=12 stations around the circumference of the pipe for 
which there are experimental values of On, 1.€., there exists a 0% 
corresponding to F(Qi ), i=1,-,N. O¢ is considered to be the experi- 
mental value if F( D;) is considered to be the smoothed value. It is 
assumed that: 

error = @; = O¢ — F (Dc) 2.4.3 
It is necessary to minimize the sum of the square of the errors in order 
to have smoothed results and experimental results within reasonable 


aguocmenth yepemetor es 


€=F(e)* 
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is to be minimized, This means: 


\ Fe 

gE = 0 
= 0 
2£= 0 

7 2.4.4 

26 = © 
g£ =O 
© 

QE- 0 
va 

3 O 


Equations 24.4 are solved in the following manner: 
Ci= Oc- F( Gz) 


= O-La+bsing+ ccosgd+dsin 26+ ecos2g + Fsin3¢ +9COS 3d] 2.4.5 
E= Zler) = Z [Oe* — 206) +[ (Be) ]*} 
a6 = 25 G SF(d + 2F(Di) SEG) = O 
3 ‘st OK K 


where K =a, b, c, d,.e,.f, g 


2.4.6 


By substituting the various values of K in equation 2.4.6 the following 


equations are obtained: 


SE = - 2200+ 22 F(DjEO 


o€ BAG 
Ia 

a€ - ~ 220% sINGe + 2E F(@e) siINGDe = 2.4.8 
95 = - 220° cOSGi + 2EF(Gi)cosGe = O 2.4.9 
Cc 

o6 = - 2EG¢ SIN 2h: + 2E F(Gi) SIN 2G =O 2.4.10 

26 = ~ 220; COS 24i + 2E F(Gi)cos 26 = O 2.4.14 
ec 
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OE = = 2ET0 SIN 2Gc+ 2EFGr)SNIG =O . .. 2.4.42 2. 
2 la 
2€ = - 22% cos 3Ge4+ 2EE(Gr) COS BGT FO 2.4.13 
9 


Shifting the first terms to the right hand side and substituting equation 2.4.la 
for F(Gi) in equations 2.4.7 through 2.4.13, the results displayed in 
matrix form are: 

(a) [B) = [c} 2holy 


where: 


Ml 589; oe, 5829; Sc2 253Q; £034; 
sf; <t8¢¢; esis <afe2g. a 28 ie g. wer 


cH «= EcPsG; Eee Schs2¢. soheag, ScPs3f; <=chc3P 


Won 


[al = kes2h ss2¥s¢- <s2 ¥ ss22g;  Es2dc ran $52ds3¢, = <286 3th 2shel5 
c2g, <c20s, ih $¢2 aa hae 2c2fs3f, <sc2Pec 
5330, 5s3s¢; 553 4830526 $83%c2¢. 8234, sf 


°36; sc3Asf, =c3 “i 5.034526; ee Sc3fs3f. =c23 





(In this matrix only, s is used for sin, c is used for cos.) 


224.16 


bo 

me 

ll 
mnMmoROAo’ & 


20% 
20i sing; 
STi cosPy 
[cj= | Sdjsin2gi. 2 ohel7 
20 ( cos2gij. 
20 { sin3di 
20 { cos3Gi. 
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When the number of stations around the circumference is even and the 
stations are equally spaced, it can readily be shown that all non-diagonal 
elements of the A matrix in equation 2.4.15 become zero, In this case the 
seven equations with seven common unknowns become seven equations, each 


with only one unknown, as follows: 
GEN) Sear 


be sintGhi = Sa SINZzGe 
@ 2 Costs 2 OEOSH Ge 
AS sint2ge = Soe SIN* 2G PEARS, 
C Sceos? 2ZGi = Sq" cos* a OO 
Pe SINeCQ! = Sov SIN 3G 
gGrecos* She 2 oecos 29 - 
By a similar derivation, the following equations are established for 
the seven coefficients used in determining the smoothed values of Tay 
[a] [B'} = [c'J 264019 


where [A] is the same as shown in equation 2.4.15, 


k 
1 
m 
[B’ | = n 20220 
) 
p 
q 
2Ua 
STi sing), 
; 2 ti.cosGi - 
| = 2; sin2gi. Lelpeew 
ST: cos2gi 
Ui sin3¢i 
Equations 2.4.17 become: 
kN = st 
lising = 2tsingi 
msicos$ = Sticod#i 
nZsine2g. = &tisin2g: 20.22 
0 Zcos<2g. = eticos2; 
p? sine3g, = etisin<3¢; 
Qs cos©3Q 2 eticos<3Gi 
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Finding the seven unknowns for Un and the seven unknowns for Txyby use of 
equations 2.4.18 and 2.4.22 will be referred to as the short method in 
the remainder of this thesis, Finding the unknowns for On and Txyby use 
of equations 2.4.14 and 2.4.19 will be referred to as the long method in 
the remainder of this thesis, 

Once the two sets of seven unknowns have been found by either the 
short or long method, they are substituted into equations 2.4.la and 
2ehe2a to determine the smoothed values of Onand lees 

THEORY FOR EQUILIBRIUM CALCULATIONS 

Equilibrium of the system was established in several steps. The 
first step was to calculate the sum of the forces and moments at each 
rosette ring cross section. The x’ direction was taken as the direction of 
flow in the pipe, the y’ direction as the wi from the center of 
the pipe toward rosette number 12, and the z‘ direction as the direction 
from the center of the pipe toward rosette number 3, as shown in figure 


Zeeks 


yf 


Figure 2.5.1 
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The second step in the equilibrium calculations is to reorient each 
individual cross section coordinate system to a standard coordinate 
system. The standard coordinate system aligns the x direction with the 
fore to aft axis of the ship, the y direction with the vertical upward 
axis of the ship, and the z direction with the port to starboard axis of 
the ship. The reorientation is obtained using actual dimensions and 
angles taken from the system blueprints and standard coordinate re- 
orientation methods, 

The third step in the equilibrium establishment consists of making 
the segments between rosette rings into free bodies and summing the 
forces and moments acting on the segment. In addition to the forces and 
moments calculated in equations 2.5.1 through 2.5.6, the forces due to 
the fluid in the pipe segment are added to obtain the equilibrium forces 
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and moments, The signs of the forces and moments at the end of the 
segment in which the fluid enters are reversed because the values cal- 
culated in equations 2.5.1 through 2.5.6 are the forces of the adjacent 
segment acting in that segment and not the forces af the incoming segment 
acting in the adjacent segment as are desired for equilibrium calculations, 
Theoretically the sums of the forces and moments on each individual 
segment should be zero. 

Only step one utilizes the computer program developed, Since each 
rosette ring is oriented differently with respect to the standard 
coordinate system chosen, more time would be rome to adjust the program 
to calculate the different reorientations for equilibrium calculations 


than to calculate the reorientations by hand. The same applies to step 


three, Therefore, steps two and three are performed by hand. 


3. Explanation of Computer Program 
SRE. ae eee «© INTRODUCTORY COMMENTS 

This section presents a full explanation of the digital computer 
program developed for evaluating the strain gage data. A listing of the 
various symbols used in the program, but not defined in the program 
explanation, is found in Appendix B, A complete listing of the program 
appears in Appendix A, The statement numbers referred to in this section 
are the six digit letters following each line of the program listing.. 

The program language employed in this program is CDC FORTRAN 60. 

All computer runs were made on the Control Data Corporation 1604 computer 
and associated equipment located at the U. S. Naval Postgraduate School 
at Monterey, California, 

The program developed is divided into 12 parts. The following sub- 
sections describe these parts. 

PART 1, DATA READ IN 

Part 1 of the program reads Aeros quantities that remain constant 
during the entire program, These quantities most likely would vary if 
used for another piping system, but would still remain constant for that 
particular systen. 

The quantities read in are N, the number of experimental depths at 
which readings were taken; E, Youngs modulus for the pipe material in 
millions of pounds per square inch; hb, Poisons Ratio for the pipe 
material used; TH(1) through TH(8), the average wall thicknesses for rings 
one through eight; RA(1) through RA(8), the pipe radius for ring one 
through ring eight4; and Z, the antic between two strain gage rosettes of 


each ring in radians. Also included in this part are the commands which 


Aohese values of radius and thickness were obtained by averaging 
the thicknesses and radii at the various measurement points for each ring. 
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initialize all AM, BM, and GM registers at the start of calculations 
for each ring. 

Statement number 000280 sndicatelithe number of rings of strain gage 
rosettes for which data is to be computed, This is indicated by the 
number of times the "DO LOOP" is used. 

PART 2, COMPUTATION OF BEST SLOPE 

Part 2 of the program reads “he actual strain gage readings in micro- 
inches per inch from data cards at the end of the program and by the 
method of least squares determines the slope of the line connecting the 
data points from 0% to 100% of test depth. The slope ‘is given by the 
value of B(M) for each Bee. The actual strain at zero percent depth is 
given by the calculated intercept, identified in ee program by AZ. The 
zero percent depth intercept is not used in this program, however, 

An indication of the accuracy of the data can be determined in this 
‘ part of the program. One of the quantities calculated in this part, 
identified in the program by DIFFS(1), gives the difference between each 
actual strain gage reading and the calculated best line fitting value at 
each reading depth. The magnitude of DIFFS(I) will indicate the accuracy 
of the data, 

in Weture applications DIFFS(I) could be used to expand the program 
to accept or reject each and every data point by inspecting the magnitude 
of DIFFS(I) compared to the values of DIFFS(I) for a particular gage at 
other depths, The program at present does not reject any individual data 
points. If one or more points were to be rejected, the program could 
be altered to compute a modified slope and intercept. 

The slope, B(M), in this program is given in strain per percent of 
test depth, 


Part 2 of the program is run three times before going on to part 3 of 
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the program in order to have all three strains available at once when 
computing the principal stresses for each rosette. 
PART 3, COMPUTATION OF PRINCIPAL STRESSES AND AXES 

Part three of Lhe program utilizes. the strain per percent of test 
depth calculated in part two of the Sropthaweee each element of the 
rosettes to determine the two principal stresses and direction of 
maximum principal stress taking the direction of flow within the pipe as 
the zero angle, The values of the two principal stresses are given in 
pounds per square inch per percent of test depth. 

Statements 000860, 000870, and 000880 check to see if any of the 
individual element strains are zero. If any Tae nloaall strain is zero, 
the program jumps to statement 001020 which sets the two principal 
stresses and angle equal to zero. 

Statements 000890 through 001000 calculate two principal stresses 


using the formlas 


Tr 13 « = E/etes + 4d VE -€3)> + [2 €2 ~it€2)J}* ] 
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and the angle using the formula 


Angle =1 are tan Zé, -(€+€s) 
a €i—-€3 

Part 2 and part 3 of the program up to statement 001050 is run twelve 
times before continuing further,. The number of times these two parts are 
run is determined by the number of rosettes per rosette ring. For the 
data evaluated, there were twelve rosettes per ring. 

Statements OO01LO60 through 001090 orient the angle of the maximm 
principal stress with respect to the direction of flow through the pipe. 
This is only necessary when one of the outside legs of the rosette is not 


oriented in the direction of flow. The rosette orientation as indicated 
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in figure 1.2.1 has caused the necessity in this particular case for 
including statements 001060 through 001090. 

A print out of values of Gmax,Gmin and Op for all twelve rosettes 

in a ring is called for in statements 001100 and 001110. 
PART 4, DECISION REGARDING MBTHOD OF SOLUTION 

Part four of the program determines whether the short or long solution 
to solve for the smoothed values of Gn and Txy mst be used, 

In subsequent sections of this thesis the use of the word inoperative 
may have one of two meanings. In one meaning,it is used to represent any 
rosette element for which data was not: recorded’, In the other meaning, 
it is used to indicate any rosette element with a zero slope as determined 
by the method of least squares, The second meaning was a decision made 
by the writer based on a close evaluation of the data without its use, 

_ When any element of a rosette is inoperative, the value of Tmax y 
calculated in part three, for that rosette is set equal to zero. This 
part of the program checks to see if any value of Omafor a ring is zero. 
If any one value of Umax is zero, the program jumps to the long method of 
solution, If none of the values of Tmax for @ ring is zero, the program 
goes to the short method of solution. 

PART 5, DETERMINATION OF FOURIER COEFFICIENTS ~ SHORT METHOD 

Part five of the program determines the values of the coefficients 
used in the determination of the smoothed values of Gn and Txy by means 
of the short method as defined in section 2 of this thesis, The computer 
space names are used in both parts five and six in determining the same 
items, 

When data has not been recorded for a rosette element, the 
strain readings entered on the program data card for that element are 


zero. This is done in order to keep proper data card sequence in the 
program, 


Statement 001230 changes the angle &from degrees to radians while 
statements 001240 through 001260 solve for the experimental values of (Gn 
and Txy. -Phi, statement number 001280,is the angle of the rosette from 
rosette number 12, Statements 001290 through 001400 compute the values 
of the sines, cosines, sines eoared and cosines squared of the single 
angle phi, twice the angle phi and three times the angle phi. The value 
of AM(1) in statement 001410 indicates the number of rosettes in each 
ring while statements 00120 through 001470 sum up the squares of the 
sines and cosines of the single, double and triple angles, phi. These 
last seven statements are the main diagonal elements of the 7 x 7 mtrix 
shovm in theory (equation 2.4.15). Statements 001480 through 001610 
calculate the values of the two 7 x 1 mtrices shown in theory (equations 
204.17 & 2.4.21). Statements 001620 through 001650 change the values of 
Gvand Txy into values needed in making the polar plots. The last number 
of the "DO LOOP" which connects statements 001220 through 001660 is 
determined by the number of rosettes in the ring. 

Statements 001690 through 001780 are name changing statements. 

Statements 001790 through 001850 solve the seven equations for the 
seven unknown coefficients for the smoothed values of Gywhile statements 
001860 through 001920 solve the seven equations for the seven unknown 
coefficients for the smoothed values of txy. These 14 equations are 
equations 2.4.18 and 2.4.22 in section two covering theory.. 

The remaining statements in this part are all print and format state- 
ments to give a print out of various items computed in this part, 

When part five is used in the program, statement 002010 is used to 
bypass part six of the program and go directly from part five to part 


SsSéevene 
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PART 6, DETERMINATION OF FOURIER COEFFICIENTS - LONG METHOD 

Part six of the program determines the values of the coefficients used 
in the determination of the smoothed values of Gn and xy by means of the 
long method as defined in section 2 of this thesis, 

Statement 002080 is used to determine if any rosette in the ring is 
inoperative. If any rosette is inoperative, statements 002090 through 
002150 are used and the program goes to the end of the "DO LOOP", This 
is used to prevent plotting an experimental value of Gy and Txy on the 
output plots when the rosette is inoperative. 

Statements 002160 through 002290 are identical to statements 001230 
through 001360 of part five. Statements 002300 through 002480 compute 
the combinations of sines and cosines of the single, double and triple 
angle phi, The items identified by AM(1) through AM(28), statements 
002490 through 002760, are the various elements of the 7 x 7 matrix shown 
in section two (equation 2.4.15). Statements 002770 through 002830 are 
the elements of the 7 x 1 matrix shown in theory (equation 2.4.17) while 
statements 002840 through 002900 are the elements of the 7 x 1 matrix 
shown in theory by equation 2.4.21, Statements 002910 through 002940 
are identical to statements 001620 through 001650 in part five of the 
program, The last number of the "DO LOOP" connecting statements 002070 
through 002950 is once more determined by the number of rosettes in each 
ring. 

Statements 002980 through 003460 identify all elements of the 7 x 7 
matrix individually (Equation 2.4.15) as can be noted by the double 
subscripts. In turn, statements 003500 through 003520 do the identical 
identification for the two 7 x 1 matrices (Equations 2.4.17 and 2.4.21). 
Statements 003530, 003540 and 003550 are the first three items needed in 
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the CALL GAUSS 2 statements, The two CALL GAUSS .2/ statenents, numbers 
re 


i 
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003560 and 003610, call the subroutine GAUSS 2 into use. GAUSS 2 
solves for the unknowns in the matrix equations and is used twice to 
solve for the unknown coefficients of both Gy and xy. 

The various other statements in this part of the program are print 
and format statements which print out various quantities computed in this 
part. 

PART 7, PREPARATION FOR GRAPHICAL OUTPUT 

Part seven of the program is directly related to the subroutine DRAW, 
It is merely a series of statements which set up the labels placed on 
figures drawn by the program. More information on the titles and the 
requirement for them can be obtained by consulting Appendix D. State- 
ments 003750 through 003900 are a computed GO TO series for numbering the 
rings appropriately as computer draws the figures, 

PART 8, COMPUTATION OF SMOOTHED VALUES OF AXIAL AND SHEARING STRESS 

Part eight of the program computes the smoothed values of Un and xy 
around the entire pipe diameter and converts the values into polar 
coordinates for plotting purposes. The 360 degrees of the circle of the 
pipe circumference is divided into 628 segments for plotting purposes. 
Since the computer plotter used for figure drawing can only draw vertical, 
horizontal and forty-five degree lines, it was necessary to divide the 
pipe circumference into as many convenient increments as possible. The 
value 628 was chosen since the computer utilizes radians, not degrees for 
angles, and 628 is a convenient number with which to divide 2Tr radians 
and 900 is the limiting number of points that can be plotted by the sub- 
routine DRAW, 

Statement 003970 divides 2TT by 628. Statement 004000 computes the 
angle, Statements 004010 through 004060 compute the sines and cosines 


of the angles needed in the formlas for Gvuand ty. Statements 004070 


ad 
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and 004080 compute the smoothed values of Uw using the formula 
Iv=A+bsinG+ccosfrdsin 26 + CCcoseGrFfSin 36 + g cos SO. 
Statements 004090 and 004100 compute the smoothed value of txyusing the 
formla 
ty = k+ Ising + mcosg+nsinZG + OCOS 2h + PSIN 3G + 9COS 3. 
Statements 004110 through 004180 transform the values of Quand txy into 
values that can be utilized in polar plotting. Statements 004110, 004120, . 
004150, and 004160 are used for drawing reference circles on the figures, 
These reference circies indicate the zero stress ring. Any reading inside 
the zero stress circle for Qycan be taken as a compression, while any 
reading outside the zero stress circle for 0y can be taken as a tension. 
For the plots of Tey, any reading inside the zero stress circle can be 
taken as a negative shear, while any reading outside the zero stress circle 
can be taken as a positive shear, A positive shear is defined as a shear 
caused by the portion downstream of the rosette ring acting in the direction 
of increasing rosette numbers on the upstream portion. A negative shear 
is defined as a shear caused by the portion downstream of a rosette ring 
acting in the direction of decreasing rosette numbers on the upstream 
portion. The number 100, chosen for plotting the zero stress circle for 
Gy and 40, chosen for plotting the zero stress circle for Tey» were chosen 
to keep all plotted values positive and to show as much clarity on the 
figures as possible, If, in future applications of this program, these 
scales are not convenient, change the numbers 100 and 40 in statements 
001620-001650 inclusive, 002910-002940 inclusive and 004110-004180 
inclusive to more appropriate values. 
; Jot seat “PART 9, GRAPHICAL OUTPUT 
Part nine of the program consists of three statements which call the 


Subroutine DRAW for plotting the experimental and smoothed values of 
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Oy around the circumference of each pipe section at a particular rosette, 
a title change statement, and three more statements which call the Sub- 
routine DRAW for similar plotting of Ty values. The CALL DRAW statements 
occur in groups of three; the first one plots a circle of zero stress, 
the second one plots the smoothed points, and the third one plots the 
experimental values, Since compression stress appears as a negative stress 
while tension stress appears as a positive stress, a polar plot could not 
be utilized unless the position of zero stress was far enough away from 
the center of the plot to prevent the stresses from being plotted 180 
degrees out of position. The proper scaling of the plots must be done 
before the program is utilized. Appendix D , the Subroutine DRAW. writeup, 
gives full details on both scaling and what should be placed in a CALL 
DRAW statement. 
PART 10, COMPUTATION OF FORCE AND MOMENT COMPONENTS 

Part ten of the program computes components of the three forces and 
the three moments at each ring which are used for the equilibrium check. 
The x-direction is taken as the direction of flow in the pipe. The y- 
direction is taken as the direction from the center of the pipe toward 
rosette number 12, The z-direction is taken as the direction from the 
center of the pipe toward rosette number three. XFOR is the force in the 
x-direction and XMOM is the moment in the x-direction. YFOR is the 
force in the y-direction and YMOM is the moment in the y-direction. 
ZFOR is the force in the z-direction and ZMOM is the moment in the z- 
direction, The three force components and three moment components are 
used in the manual computation of the three equilibrium forces and three 
equilibrium moments. 

thy PART 11, GAUSS 2, SOIUTION OF SIMULTANEOUS EQUATIONS 


Part eleven of the program is the Subroutine GAUSS 2. This is a 


standard subroutine in the computer subroutine library at the U. S. Naval 
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Postgraduate School. The computer facility explanation of this sub- 
routine is Appendix C of this thesis, 
en~)» PART 12, DRAW — GENERAL PROGRAM FOR GRAPHICAL OUTPUT 

ns twelve of the program is the Subroutine DRAW, This is also a 
standard subroutine in the computer facility subroutine library at the 
U. S. Naval Postgraduate School, The computer facility explanation of 


this subroutine is Appendix D of this thesis, 
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4. Data Evaluation 
GAGE ORITENTATION 

The gage orientation is one of the important items to which careful 
thought should be given in order to obtain the most meaningful data. The 
most important consideration concerning gage orientation is the number of 
rosettes per ring and the spacing between rosettes. As previously discussed 
in section 2, it is easier to obtain smoothed values of ON and ‘Ty when 
an even number of rosettes equally spaced are employed because equations 
with only one unknown are more readily solved than several equations with 
several unknowns. The actual number of rosettes per ring employed in the 
SSBN 640 instrumentation appears to be a good number for the size piping 
system to be analyzed. On future instrumentation projects the size of 
the piping system should be considered in determining the even number of 
gage rosettes to be used. 

In the system analyzed the rosette orientation was shifted 45° 
alternately between even and odd numbered gages (See figure 1.2.1). 

This shift in orientation caused minor difficulties in obtaining the 
proper angle between the maximum principal stress and direction of flow 
which had to be checked very closely; it is also a source of possible 
misunderstamiing. The identical results could have been obtained with all 
rosettes oriented identically and thereby eliminating these minor diffi- 
culties, lt is highly foCoulichTee for future instrumentations that all 
rosettes be oriented identically. 

It might be beneficial to consider the use of T- DELTA rosettes for 
future instrumentation projects, In several rosettes, one leg indicated 
erroneous readings. If a fourth leg of the rosette had been available, 
the rosette might have produced meaningful results not otherwise obtain- 


able @ 
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ACCURACY OF READINGS - DOCK TRIAL DATA 

The readings taken during dockside hydrostatic tests were recorded 
to the nearest 10 Ain/in. Since it is possible to read an SR-4, Indicator, 
Type N, more accurately than the nearest 10 }in/in, the accuracy of all 
readings is not considered to be the best possible, I1t is desirable in 
future tests where an SH-4 Indicator,Type N,is used to read or estimate 
strain gage indications to l fin/in. Seven pressure increments were 
recorded during the tests. 

One of the items printed out in Part 2 of the digital computer 
program is DIFFS(I), the difference between each actual strain gage 
reading and the least squares calculated best line fitting valup j In 
evaluating the actual data, any value of DIFFS(I) greater than 10 Min/in 
is considered to indicate an inaccurate reading. With this as a basis 


of inaccuracy, the following statistical tabulation of the dat has been 


made, 

Table 1 
Total number of readings taken 1930 
Total number of readings with DIFFS(I) greater than 10 sin/in 19 
% of total in error 2.52% 
Total number of gage elements 288 
Number of inoperative elements Del. 


Number of elements with one reading more than 10 win/in in error 15 
Number of elements with two readings more than 10 w&in/in in error 9 
Number of elements with three readings more than 10 fiin/in in error 4 
Number of elements with four readings more than 10 A&lin/in in error_1l 


Total number of inaccurate or inoperative elements 1,0 
% of total elements that were inaccurate or inop@rative 13.89% 
Total number of rosettes 96 
Number of rosettes with at least one inaccurate reading or 
an inoperative element Ring 1 4 
Ring 2 3 
Ring 3 L, 
Ring 4 S| 
Ring 5 3 
Ring 6 es 
Ring 7 2 
Ring 8 Zs 
Total 30 
% of total rosettes inaccurate or inoperative 31.25% 
Total numper of rosettes inoperative due to faulty element 7 
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~ of total rosettes inoperative due to faulty element 9.375% 
Ring 6 had one rosette with all three elements faulty 


The statistical percentages shown in Table 1] have been interpreted 
in the following manner: 

a, Faulty elements caused 9.375% of the rosettes to be unusable 
in evaluating the systems, Although this percentage is considered to be 
high in obtaining good experimental results, it is not excessive in 
evaluating this system. 

b. The percentage of inaccurate rosettes caused by a reading 
error of 10 Auin/in, 21.875%, does not have any real meaning because 
only one, two, or even three inaccurate readings might not have caused 
the slope of the strain depth curve to be enough in error to make the 
rosette*give a faulty overall stress calculation. 

ec. The percentage of inaccurate readings, 2.524, is well within 
the experimental accuracy expected, if noteven less than could be expected 
from a test performed under the conditions of this test. 

PRINCIPAL STRESS RESULTS - DOCKSIDE HYDROSTATIC TESTS 

The Bureau of Ships Shipbuilding Specification,MILT 16420G, states 
that the maximum allowable stress in 70-30 copper nickel piping under 
cyclic loads is 15000 psi. In publications on material properties, the 
yield strength of 70-30 copper nickel is 25,000 psi. The maximum value 
of Oma, at any rosette at 100% test depth calculated from the dockside 
hydrostatic test data was found to be slightly less than 10,000 psi. 
This indicates that the maximum principal stress does not reach the allow- 
able stress and is well below the yield strength at any instrumented point. 
Since the rosette ring locations were chosen to be representative of the 
entire system, this is sufficient proof that in all probability the 
stresses within the piping system do not exceed the allowable stress 


limit of the material, It is also proof that the stress levels were 
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‘ 
sufficiently low to assure elastic action, as is assumed in the analysis 
of the data, 
COMPARISON OF EXPERIMENTAL DATA TO SMOOTHED CURVES 
The sixteen figures appearing in Appendix E are polar plots of 
smoothed axial and shearing stress around the circumference of the piping 
at each rosette ring and the experimental values of the axial and shearing 
stress at the individual rosette locations. The circle ten units in 
radius on the Gy plots and four units in radius on Txy plots is the circle 
of zero stress on these plots. For the Gn plots stresses inside the zero 
stress circle are compression while those outside the circles sre- 
tension,. For the 1e plots stress inside the zero stress circle is 
shearing stress in the negative direction while stress outside the zero 
stress circle is positive shearing stress. The scales on the plots are 
arbitrary to remain within the unclassified security classification as 
noted in section 1. The stresses shown are plotted per percent of test. 
depth, a | 
It is readily apparent by close examination of the sixteen figures 

that the experimental values are in close agreement with the smoothed 
curves. (The experimental points are marked by an x while the smoothed 
values appear as a closed curve.) The x at the center of some plots:is 
to be ignored as it is only a rosea, step in having the CDC 1604 
computer draw the plots. It will also be noticed that some plots do not 
have all twelve experimental points shown, The missing experimental 
points are an indication of the location of faulty rosettes (inoperative 
elements) as no data from fully operative rosettes was rejected in this 
evaluation. Had the missing rosettes been operative, the shape of the 


curves might have been altered. 
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EQUILIBRIUM RESULTS FROM DOCK-TRIAL DATA 
Equilibrium calculations have been performed on the segments of piping 
between rosette rings l and 2, 3 and 4, and 7 and 8 Calculations were 
not performed on segments connecting any other two rings because unknown 
forces due to foundations, pipe hangers, or other constraints exist 


between these segments. Table 2 gives a tabulation on the three segments 


analyzed. 
Table 2 : 
Fx|ry [Pz [ux [My melee ee eel a 
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(Resultant moments are about lower numbered end of segments ) 





Theoretically, for equilibrium to be satisfied, all values of Frx, 
Frys Fess Mes Mey, and: Mezshould be zero. Although only one of the 
eighteen resultant values is zero, the other values are not excessively 
high. Any one of the following influences could have..caused the non-zero 
values. 

a, The distances between rosette rings that were used in 
calculating the resultant moments were distances measured from the piping 
blueprints using the approximate locations of the rosette rings. Since 
exact values of these distances were not available, these were the best 
values available to use even though they could be a source of error in the 
calculations, 

b. The values of pipe radius and thickness used in calculating 
the forces and moments at each rosette ring were average values calculated 


by averaging the values of thickness and radius at each rosette within the 
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ring. In addition the values of radius at each rosette reported in 
Propulsion Plant Test Form 721-5 were obtained by measuring pipe diameters 
and dividing by two, which is not a true value of the radius considering 
that the circumference of the pipe is not a true circle, Since outside 
diameter measurements for rosette ring number 1 were not recorded, the 
value of radius of ring 1 used in force and moment calculations was 

the average of the other seven rings. 

c. The-value of the forces exerted by the fluid within the pipe 
on the inside surface of the pipe was determined using average values of 
radius and wall thickness. Using average values instead of exact values 
could have once more caused errors. 

d. Errors in the Fourier coefficients caused by errors in strain 
gage readings could also have contributed to errors in the forces and 
moments calculated at each rosette ring. 

Considering all the possible sources of error in the equilibrium 
calculations, it is reasonable to expect that all the resultant forces and 
moments would not be zero. However, some of the values in table 2 are 


much larger than could be considered reasonable for the method used, 
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5. Conclusion 
a ae VALIDITY OF RESULTS 

As indicated by the results shown in this thesis, the method 
developed can readily be used to obtain stress information concerning 
the system. However, the validity of the equilibrium results will depend 
largely on the accuracy of the strain gage data and other data needed for 
the equilibrium calculations, 

In order to increase the accuracy of the data and thereby increase 
the validity of the results, the following recommendations are submitted: 

a. Appoint a project engineer with experience in piping theory 
and strain gage measurements to have overall control of the instrumentation 
of the system, data taking during all tests, and reporting of results 
obtained from the data. This one person would supervise the project from 
the start of instrumentation to the submission of the final report. Having 
one man in charge of the project from start to finish will greatly reduce 
the chances of faulty communication between people working on the project, 
reduce the chances of necessary data items being omitted from the data 
taken, and reduce the chances for error in the actual data taken. 

b. Utilize automatic scanning and recording equipment to record all 
strain gage data during actual tests. Automatic recording equipment can 
observe strain gage reading much more accurately than equipment read 
manually and in much less time. In addition, human recording errors could 
be eliminated. A rapid record of strain gage reading at any depth could 
also be taken using the automatic scanning and recording equipment. 

c. An important aid to the person analyzing the experimental 
data from these tests would be notes taken during the period that data 
is being taken. Some of the information that would be useful as notes 


is any unusual conditions that occurred during data taking, comments 
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concerning any data that might give an immediate indication as to being 
inaccurate or unreliable, the conditions that existed during the periods 
between the taking of data,, @tCee The project engineer should be 
observing all tests to take his own notes as well as having notes taken 
by all technicians involved in the tests, 

d. Additional data to be taken during the tests should include 
actual measured distances between two adjacent rosette rings for use in 
more accurate equilibrium calculations and the location of the welded 
seams for each section of pipe. 

By further refinement of data taking methods and actual data taken, 
the method of evaluation developed in this thesis can be used as a tool 
in analyzing any piping system instrumented in a manner similiar to the 
Main Seawater Piping System on SSBN 640. 

Certain refinements to the present method of evaluation have been 
considered; however, sufficient time is not ,available to permit these 
ideas to be included in this thesis, The basic content of these ideas 
is listed in the following paragraphs. 

Part 2 of the program could be expanded to include a method which 
would investigate the value of DIFFS(I). This would be a two stage 
process, The zero percent depth intercept, a; slope, b; and values of 
DIFFS(I) would be found as they are in the present program. If DIFFS(TI) 
were greater than a predetermined limit value, the data which produced 
that DIFFS(I) value would be rejected and new values of zero percent depth 
intercept and slope would be found, An addition of the type proposed 
would improve the accuracy of all resultant calculations. 

The present theory might be further developed to find a method to 


make some use of strain gage rosettes which have one or two inoperative 
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elements, When the center element is the only one inoperative, some useful 
informtion may obviously be obtained from the rosette. If the orientation 
is that of the even numbered rosettes as shown in figure 1.2.1, elements 
A and C operative and B inoperative would permit evaluating the axial 
stress, If the orientation is that of the odd numbered rosettes as shown 
in figure 1.2.1, elements B and C operative and A inoperative would permit 
evaluating a shear stress. Beyond this, if rosette failure were other 
than as indicated above, there is probably some statistical way of making 
use of the data afforded by the operative elements. However, the theory 
of such usage remains to be worked out. 

These two refinements added to the analysis method might produce more 


meaningful results from the dockside hydrostatic test data, 
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244 ALLE ST Le SLGXXL12 be TAU Y (ee Tee ee (2 ee ( ES 
SOON 2 yee cee a NO! EZ TS TSE VERY ee yen Ley ee Se ee | 6 ence ee 
QUT 2) 3 Ponte la 37 1s VSR ( 6 30 pe XSN7 1639) 9% oN SA ls X 1 O52 oy 6 ee 
BT ((639 )0¥ Woi/(635°:) el (Sls TAUS pe MEOR(S js yeeeMiS ) 6 oe yk eee 


[een Sant or PAST Te PAR” < KEAGS 


oy Case orl (7) 
] Thy DATHe 


rEAD 33959 Nek ot 
399 FORMAT (L14s2F 1545) 
~H(1)=92e625 
rP4(1)2=664C 
TH( 23292634 
RA(2)=6e394 
TH(3)=0 2629 
RA(3)=6e409 
TH (4)=02633 
RA(4)=6e374 
TH(5)=0 2428 
RA(5)=6e416 
T4H(6)=62640 
RA(6)=6e354 
TH(7)=0 0423 
RA(7)=6e411 
TH(8)=02423 
RA(8)=6e410 
DO 6&5 K=198 
pO 3 JT=1928 
AM(1)=0. 
DO 4 T=1ls7 
TOL =O 6 
4 CM(1)=O06 
2=32e14159/6.6 
THIS IS THE END OF PART le 
14 60°15 E=I412 
THIS IS THE START OF PART 26 
PROGRAM [TO DETERMINE THE STRAIN PER PERCENT OF DEPT FOR EACH 
INDIVIDUAL GAGE. 
57 (po a6 M=— 1a 3 
R=AD 1009 (XCI1)sYCT)sT=19N) 
1G. FORMAT(14F5.9) 
CO 996 IT=1 5N 
XSAVECTIEACT) 
991 yY SEVEC 1 1 YA 1) 
Voie = Oe 
SUumMX = Ge 


Les) 


SUMY = Oe 
Skee se (Che 
VW)WAXX = 6 
Mile 2 4 


Le 3 1 = ds. i 
SINAX=SUMX4+X( 1) 
Sy = SUMY+Y¥ C1) 
Sit Aw Cs AX Rae i A eek a 
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PART 2 READS THE DATA AT THE END OF THE 


009370 
000 380 
000390 
900400 
000410 
000420 
000430 
000440 
000450 
000466 
000479 
900480 
9nNN4&9c 
000500 
000510 
000520 
000530 
009540 
C00550 


NAN 


37 SUMXYVYESUMXY4 (X(T) #YC71)) 
DFNOM = XN¥SUMXX=-SUMX#SUMX 
ANUM = SUMY &® SUIMXX -— SUMXK *® SIJYXY 
AZ = ANUM / DENOM 
RNUM= XN#SUIMXY=$SHUMXKESUMY 
B(M) = BNUM / DENO™ 
DO 40 J = 1s5N 
YC = AZ +B(M) *xX(J) 
DIifpe= YCJ) =)YC 
4C YSUM = YSUM + DIFFROILFF 
S2 = YSUM / (XN —- le) 
SI1G=S2/DENOM 
SRA=SORTF (SUMXX*SIG) 
SRB=SQRTFU(XN*¥SIG) 
PRINT 9975 AZs RB(M)s SRA» SRB 
997 FORMAT(24HCOC STRAIGHT LINE FITTING// 25H USING THE FORM Y=BX + As 
1/ 4H A = F1l4.8/ 7TH R(M) = F14,8/ 26H STANDARD DFVIATION OF A = E14 
208/ 29H STANDARD DEVIATION OF RBI(M) = £14¢8) 
801 SUMFIT=.0 
DO 90 I=lsN 
YFIT(I) = BIM)*X(1)+AZ 
DIFFS(I)HVYFITC(II-YCI) 
90 SUMFIT=SUMFIT+ DIFFS(I) 
995 PRINT OS sSUMFITs(X( 1) sYSAVE(T) sY(C I) sVYFIT(I) sOILFFS(I) sI=1leN) 
99 FORMAT(SHOSUMFIT =E12¢5///10X1HX1L4X1HYSX1LOHY TRANSPOSEBX4HYFIT11X5H 
HDTEESS 771 5E 1565) ) 
58 CONTINUE 
THIS IS THE END OF PART 2 AND THE START OF PART 3-6 THIS PART 
DETERMINES THE PRINCIPAL STRESSES AND THE ANGLE FROM THE PIPE 
CENTERLINE OF THE MAXIMUM PRINCIPAL STRESS FOR EACH ROSETTE. 
JF(B8(1)) 461941054091 
401 IF(3(2)) 402594199402 
402 IF(R(3)) 40394109403 
403 AA=E/2-6 


SA=B(1) + B(3) 
GA=601)> = 8(3) 
D=1.+U 
F=1.-U 
G=2-*B(2) - BA 
H=BA/F 


P=1l1e/D*SQRTF(CARCA+G4G ) 
SMAX (LL) =AA*(H+P ) 
SMIN(L)=AA® CHP) 
ANG(L)=005* (ATANF(G/CA )4(3014159/20) * (1 e-(CA/ABSF (CA) )))*¥1800/3014 
1159 
GO TO 15 
410 SMAX(L)=0O64 
ANG(L)=Ce 
> CONTINUE 
DO 115 L=alelles2 
115 THETA(L)=ANG(L)-45.6 
DO 116 L=2s129s2 
116 THETA(L)=ANG(L) 
PRINT So(SMAX(L)sSMIN(L) sTHETAI(L) sL=isi2) 
5 FORMAT (///6X THSMAX(LIBXTHSMIN(L) BXBSHTHETA(L)//7(3E150e5)) 
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000560 
000570 
000580 
000590 
009600 
000610 
000620 
000630 
000640 
000650 
000660 
000670 
000680 
000690 
000700 
000710 
0007290 
000730 
000740 
000750 
000760 
000770 
000780 
000790 
000800 
000810 
000820 
000830 
000840 
000850 
000860 
000870 
000880 
009890 
000900 
000910 
000920 
900930 
000940 
000950 
000960 
000970 
000980 
000999 
001000 
001010 
001020 
001030 
001040 
001050 
901960 
001070 
901080 
001090 
001100 
001110 


a ta 


AY OX 


Tle Soe idea 
WHETHER THE § 
MOST BE USED. 
DO 429 L=le1l2 
IF (SMAX(L)) 42094379420 


Ines slaxl Gs 
UT TP ohraner sc ae 


DART be 


NO CF PART 3 AND 
ACRT OR G ‘s Tae “Ses > sea 


GG 


Cr ae 
ee 


42C CONTINGE 


THIS 15 THE END CF PAPT 4&4 AND THF START OF PART Se PART 5 DETERMINES 
THE VALUES OF THE COEFFICIENTS USED IN THE THEORETICAL DETERMINATION 
OF NORMAL STRESS AND SHEAR STRESS BY THE SHORT OR EVEN NUMBER OF 
G4GFS EQUALLY SPACED METHOD. 

DO 16 L=1ls12 

TRETACL)=THETA(L) #3214159/18C, 

SIGXX(LI=ESMAX(L)I*¥COSF (CTHETA(L) ) ¥*COSF (THETA(L) ) +Si4IN(L)*#SINF( THETA 

1L))*SINF(THETA(L)) 

TASXY(L)=(SMAX(LI-SIMIN(L) )¥O00S¥SINF(2e*THETAI(L) ) 

Ku 

PHI(L)=XL*¥Z 

SI=SINFC¢(PHI(L)) 

CI=COSF(PHI(L)) 

SI2=SINF(2e*PHI(L) ) 

CI2=COSF(2e¢*PHI(L)) 

SI3=SINF(3.*PHI(L)) 

CI3=COSF(32e#PHI(L)) 

SISQ=SI*SI 

CISQ=CI*CI 

$25Q=S12*S12 

C2SQ=C12A"C1 2 

$3S5Q=S13*S13 

C3SQ=C13*CI13 

AM(1)=12.6 

AM(6)=AM(6)+SI1SQ 

AM(10)=AM(10)4+CISQ 

AM(13)=AM(13)+S2SQ 

AM(15)=AM(15)4C2SQ 

AM(26) =AM(26)+S35SQ 

AM(28)=AM(28)+C3SQ 

BM(1)=BM(1)4+SIGXX(L) 

BM(2)=BM(2)+SI*SIGXX(L) 

BM(3)=BM(3)4+CI*SIGXX(L) 

BM(4)=BM(4)4+S1I2*SIGXX(L) 

BM(5)=BM(5)4+CI2*SIGXX(L) 

BM(6)=BM(6)+S13*SIGXX(L) 

BM(7)=BM(7)+CI3*SIGXX(L) 

CM(1)=CM(1)+TAUXY(L) 

CM(2)=CM(2)+SI*TAUXY(L) 

CM(3)=CM(3)4+CI*¥TAUXY(L) 

CM(4)=CM(4)4+SI12*®TAUXY(L) 

CM(5)=CM(5)+CI2*TAUXY(L) 

CM(6)=CM(6)+S13*TAUXY(L) 

CM(7)=CM(7)4+CI3*TAUXY(L) 

XSXX(L)=(1000+SIGXX(L) )*CI 

YSXX(L)=(1000e+S1GXX(L) )*SI1 

XTXY(L) =(40 e+ TAUXY(L)) ¥CI 

YTXY(L) =(40 e+TAUXY(L) )*SI 
16 CONT INUF 

PRINT 7se(SIGXX(L)»TAUXY(L) sL=1912) 
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COL 30 
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Gorse 
901160 
SORTS 
O0 tec 
001190 
901200 
JO12iS 
001229 
001230 
001240 
001250 
001260 
001270 
901280 
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001300 
001310 
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001330 
001340 
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001360 
001370 
001380 
001390 
001400 
001410 
001420 
001430 
001440 
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001460 
001470 
001480 
001490 
001500 
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001520 
001530 
001540 
001550 
001560 
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001580 
001590 
001600 
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001620 
001630 
001640 
001650 
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001670 
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KE QRIAAT( S/S /6XY5HSO IT GYY ET TYG R TA IVY Pept 2 
le. yan (i) 

A(2se)=4%(6) 

Aldae3) =AY%i Iv) 

Al4o4)=h (42) 

A(Ss5)=A4(15) 

A(696)=A%"( 25) 

meeiesnT j=, 28 ) 

Pec ) he bene 
Bie) =38MC05 
CN(]91)=CMCi 
ZX(191)=BNC6 193 
meee yl y=oNtC ec 6] 
Press} p=: 39} 
ZX(491)=BN6491)/4( 494) 

Pree: ) = GN 541)/A1 555) 

ZXlosi)d=BN(G01)/Alo0s6) 

Peeve P-tNI ret) /AC Ts 7) 

eee PHN Gi SACL» 3) 

Pees = (nics l)/A(2Z52) 

Reece) —-ONC391)/A( 3s 3) 

XYOGol)I=HCNC eel) SAl 44) 

eee P= Hel) /Al 5s 5) 

Peeves) )=CN(Ho1)/A(696) 

meee eo. )=CNC 791) /A0 157) 

ere Gee XC ie] j)sZX'201)eZX03o1) 028145 192X091) 524055 leet 
moktiatis7/1iH ZXtlei) = EL4e¢8B//11H ZX(251) = ELl4.677 11a Zhe 
PELGeB//IIH 2X(491) = EL4GeS6//1LIH ZXt501) = EL4Sec//ilH ZXl6si) = 
Pees lr 2X isi} = €14¢8//) 

Prin LoOeK¥ (Lele 2¥ 0291) eKVY iGo ll} oXV (Gok pa XV U5 sLpsA YT lGsl i sAYl (oa 
Berlin ty7* Li XYCieL) Ei¢ed//1LIH X¥leel) = EBlSe 6/73 1 eo ee 
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1)/4( 397) 


Fils 


Metteos7tia XY¥l(oeel) = E14e8//1}H XY¥(501) = £14e8//11H XY¥(65h) = E1s 
Peorvs7ibh AY(7sl1) = El4SeBs/) 

SO 10 456 

TS §% FAP FNY OF PART S AND THF START OF PART Ge FART & DETERMINES 
feeseuesS OF THE COEFFICIENTS USED IN THE THECKETICAL DedexMiNnNAl iON 
mee vem S OF THE COEFFICIENTS USFU IN THE THEORETI GAe 
TERMINATION GE NORMAL ST&t SS AND SHEAR STRESS RY THE LONG OR 
Bea ciy SPACED GAGES AZTHOJ. 

DG 46 LFlisi2d 

TF CSMAXIL)) 449942 9449 

SLOAX(L }= 06 

TAUXY (L)=%. 

XSXX(L)=O64 

YSXX(L) FG e 

XTXY(L)=06 

WIXY(Lj}=056 

39 TO 46 

Dimer Y= THETA(L)< 36141597180, 


Stak = SAX (LY SCOUSF (THETACE) )FCOSF UTM TATL PSA NNT toe 
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TAUXY(LI=CSMAX(LI-—FAIN(L) FDC SKESENE C2 e THC TA CLD 

XL=t 

Bei (LC Y=ALEZ 

mre Siir (PHItL)) 

Pe Goorarnmit(l)) 
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002040 
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202070 
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002160 
NN02170 
2O2 180 
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5 OF OY SY me pat 
bt met mt mw ty pd tp) US ee ee FR RO ND? PB 
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GZCS=C 
$3502514*813 
Bec 3551 32C)3 
C42SQ0=C14*¢13 
AM(]TI=AM(1) 41. 
AM(2)=AM(2)+4SI 
AM (3)=AM(3)+CI 
AM(4)=AM(4)4S12 
AM(5)=AM(5)+CI2 
AM(6)=AM(6)4SI15S9 
AM(7V=AM(7)45C 
\M(8)=AM(8)4S5S1S2 
AM(9)=AM(9)4+51C2 
Afé(i10)=AM(10)4C1SQ 
AM(11)5AM(11)+4+CI1S2 
AM(12)=AM(12)4+CIC? 
AM(43)=4M(13)+S52S0 
AM(14)=AM(14)+4+S2C2 
ae }=AM(15)4+C2SQ 
M(16)=4%(16)+S13 
WEB Go ee 
AM(16)=AM(18)+S1S3 
AM(19)=AM(19)451C3 
AM(20) =AM(20)+C1S3 
AM(21)=AM(21)+CIC3 
Mt ZZ2=AM(2Z2)4+5253 
AM(23)=AM(23)+S2C3 
AM(24)=4M(24)4+C253 
BiA(25)5AM(75)4C2C3 
AM( 26) =AM(26)+5359 
AM(27)=AM(27)4+53C3 
AM(28)=AM(28)+C35Q0 


— 

iw) 
bx4 

ry 


‘4 


eR Se eee at 
M(2)=BM(02)+SIT*SIGXX(L) 
oN Siren ete IGXX¢L) 


py toy Tid) fy ey 
CG (SS (IES) De) 6) 
No Boa KR: RERQ RY fh hd NO RQ 
ys bab cw Ro KR) DN BRE ND IND 
~ CAIN 0 (5) 1D fo) Sl teh Su Je 
BY KD ey eo) (OT eT ea ag) 


Oro we > 


092409 
0C2410 
002420 
002430 
602440 
002450 
002460 
002470 
GO2489 
002499 
602500 
002519 
002520 
0025388 
002540 
002556 
002560 
C02570 
002580 
002590 
002600 
002610 
002620 
002639 
002640 
002650 
002660 
002670 
002680 
002690 
002700 
002718 
902720 
002730 
002740 
992750 
002760 
002770 
002780 
002790 
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(ead 


Wee pet (Ce eS Tee GFX OL 
Wee lor Coe Leaks porn li) 
RM (AYSEN(AP FECT AETIC“NX(L) 
AM UT \ =O" (7) tC 7 AeSIGXX(L) 
CS (lee Pyare xy (LL) 

eo (Zhe (2) +518 oa iXY(L) 

NE GT) gl ee 

CAA =e eA IXY (CL) 
eo ae iy ec ee a IXY CL) 
TMCOVECM(O) 4ST 38 TAUXY(L) 
Ct lame Pepe art aAUXY OL) 


KS KO) = ee SIT GXK(L) ) FCI 
Vee et AK CL) RSI 
XTXY(LIHC4O + TAUIYV(L)) ACI 
YI GeyHa C4 ott ATINY(L) 1 #5] 


ea) EE) 


PRINT 77s (SIGXX(L) sTALXY(L) st=)512) 


BORMagpees Om OX XI] ASHTAUXY ///(2E) 505) ) 


RCV eas 
A(t 92) =AM%(2) 
A(isS)=AM(3) 
QA(ie4)=AM(4) 
AlisS)=AM(9) 
A(is6)=AM(i6) 
ACls @)=AM(17) 
A(291)=4™M(2) 
A(2e92)=AM(6) 
A{2923)=AM(7) 
A(Z294)=AM(39) 
A(295)=AM(9) 
A(296)=AM(i8) 
N(297)=AM(19) 
A(391)=AM(3) 
A(392)=AM(7) 
A{3s3)=AM(i9) 
A(3s4)=AM(1i1) 
A(395)=AM(i2) 
A(396)=AM(20) 
Meese pt AMt?l) 
A(49i)=AM(4) 
A(492)=AM(B) 
A(493)=AM(11) 
A(494)=AM(13) 
A(4,5)=AM(14) 
Al(4s8)=AM{(22) 
A(457}=AM( 23) 
A(591)=AM(5) 
A(592)=AM(9) 
A(595)=AM(1i2) 
A(594)=AM(i4) 
A(595)=AM(15) 
A(596)=AM(24) 
A(5o7?)=AM( 29) 
4(691)=AM(16) 
A(692)=AM(18) 
A(6592)=AM(20) 
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nueoog 
302610 
MO ermes 
702838 
WIVZRSN 
062859 
2092 869 
CO2ga 
302880 
OCzere 
Coge700 
202910 
Oe Feo 
092930 
OO02940 
002950 
C2260 
9902975 
992989 
JOZ250 
903000 
003010 
003020 
003030 
003040 
003050 
003060 
003070 
003080 
003090 
903190 
02:11:60 
003120 
003130 
003140 
003150 
903160 
003170 
003180 
9003190 
003200 
003210 
003220 
003230 
003240 
003250 
003260 
003270 
003280 
003290 
003300 
0034 16 
003320 
003330 
003340 
903350 


A(694) =A" (22) 
A(695)=44(24) 
A(596)=AM(26) 
A(69r)=AM(27) 
£(792)244(727) 
A( fee) =A(1s) 
A(7s93)=4%(21) 
AC’ o4)=A%(23) 
A( 755) =4M(25) 
A(7s6)=A4A(27) 
Al 79 7)=AM(28) 
Pale 29Ss((AlJs TVs lals7) sJ=1<7) 
299 FORMAT(IS//2X5HA( Ls LISGXOHAT s 2) 4% OHAL1 93) 4K OHA 1 eG) GX6RAL 195) 4X6HA( 
LIs6)4X6HACIT so fIS// 0 7ZIV eH) ) 
DG? l=lei 
INC Ts1)=R84(7) 
Chi lsl)=CMCT) 
NN=7 
MM=] 
EP=l1eF-9 
CALL GAUSSZINNeMMo=P 9A sp BNo ZX KER) 
PRINT 7996ZX(1$1)9ZX%(291)92X(391)9Z¥ (491) 592X155] ) 92% (65 je Zeer 
7927 PORMAT(I//711H ZX(isl) = ElSe¢8//11H ZXl291) = Es4e8//711H ZAC 3s ie 
LEG eB//I1H 2X(491) = ElL4GeB//11H 2ZX(591) = ELS e38//i11H Z2ZXl691i1) = E14 
Peo/FilH ZXt7s1) = E14.8/7) 
CALL GAUSS2 (NNoit4s EP oA oCNo XY OKER) 
PRINT 75SoXYC Lol oXVOCZoldoXV (Zool oXVl(4o1) oXV 6591) oXV(691) oXY¥C 7% 91) 
Boo FORMATU/7/11IH XY(1ls1l) = El4e8//11H XY¥(2s1) = E1408//110 X17 0se eee 
TEL46¢87/11H XY¥(64901) = El14e¢8//11H XY(591) = E1LSe¢8//11H XY(65a) =]— ee 
Ceo Peek y Cis 1) = El46877) 
THIS IS THE END OF PART 6 AND THE START OF PART Te 
THE TITLES THAT ARE PUT ON THE PLOTS DRAWN BY THE 
S56 (D0 Lu? [=1s5 12 
PO? TTITEEC! p=8H 
ITITLE(1)=8H KOCHs J 
ITITLE(2)=8He Wes ix 
hit Weel 7 yack SIG VS6 
ITITLCE(S)=6H PH] 
ITITLE(9)=8HRING NO. 
GO Oo Is 3972 9303 9 304s 305 630653075 398) 55 
Pio eeeroy=3H ONE 
GeO 31> 
IPTRTEECIO =a 
GO 315 
ITITLE( 10)=8H 
GOTO. 315 
Pye 10) =sH 
GO TO 315 
ieee 6H 
GO” 1G 3125 
et e2e0 }=38H 
COMO 451-5 
Pee es i=sH 
GO 10 315 
JTITLE(10)=8H 
1 a al a Se 6 | 


PART 7 READS IN 
PROGRAM. 


TWO 
THREE 
FOUR 
ELVE 
306 Sx 
307 SEVEN 


ErGHT 
3/7/66 


396 
Bplso 
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03360 
003379 
003369 
903330 
993499 
003410 
903429 
903430 
O93449 
9063450 
993460 
003470 
003489 
003490 
003590 
0Osar¢ 
003529 
003520 
903549 
003559 
S03 S60) 
003370 
003580 
003590 
003600 
003610 
003620 
093630 
903640 
0036590 
003660 
003679 
003680 
003690 
003700 
003710 
003720 
003730 
003740 
0037350 
003760 
003770 
003780 
003790 
003800 
003810 
003820 
003830 
003840 
003850 
003860 
003870 
003880 
003890 
003900 
003910 
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‘S 
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NAD 


NAN 


C 
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THIS IS THE END OF PART 7 AND THE START OF PART 8e 


LB=4H 
LD=8H 
PART 8 COMPUTES 


THE THEORETICAL VALUES OF NORMAL STRESS AND SHEAR STRESS AROUND THE 
PIPE DIAMETER AND TRANSFORFS THEM INTO VALUES THAT ARE TO BE PLOTTED. 


26 


THIS IS THE END OF PART 8 AND THE START OF PART Qe 


Z22Z2=6.28318/628.- 

DO 26 1=1s628 

Z1=! 

PHI1(1)=Z1*22ZZ 
SO=SINF(PHI1(1)) 
CO=COSF (PHI1(1)) 
SO2Z=SINF(2e*PHI1(1)) 
CO2=COSF(2e*PHI1(1)) 
SO3=SINF(3e*PHI1(1)) 
CO3=COSF(3e*PHI1(1)) 
SIGN7( 1) =ZX(191)4ZX0291)*S04+Z2X(391)*CO+ZX( 491) *¥SO24+ZX( 591) *CO2+ZX ( 


1601)*SO03+ZX(791)*CO3 


TAUNTC TI) =XY (C1 91)4+XY0291)*SO+XY (391) *¥CO+XY (491) *SO024+XV (591) *CO2+XY ( 


1691)*SO3+XY(791)*CO3 


XSNR(I)=1002*CO 

YSNR(1)=1002¢*SO 

XSN7(1)=(1000e+SIGN7(1))*CO 

YSN7(1)=(100¢+SIGN7(1))*SO 

XTNR(1)=402*CO 

YTNR(1T) =406*SO 

XTN7(01)=(40¢+TAUN7(1))*¥CO 

YTN7(1)=(406+TAUN7(1))*SO 

CONTINUE 

PART 9 CONTAINS 


THE STATEMENTS WHICH CALL THE SUBROUTINE TO PLOT THE STRESS 
DISTRIBUTION PATTERNe 


THIS IS THE END OF PART 9 AND THE START OF PART 106 


CALL DRAW(628 »XSNReYSNRol 90sLDeoITITLE 29506950 0949492929898s09LAST) 
CALL DRAW( 628 oXSN7sYSN792 909 LBs ITITLE 9 500 950094249292 9898 909LAST) 
CALL DRAW(129XSXXsYSXXs93e9l]1 sLDeI TITLE 250 6 9504249492929898s09LAST ) 
ITITLE(7)=8H TAU VSe 

CALL DRAW(628 »sXTNReYTNRol s0eLDsITITLE 9200 29200 949492929898 s09LAST ) 
CALL DRAW( 628 oXTN7sYTN79220sLBoITITLE 9200920 0949492929898 909LAST ) 
CALL DRAW(12 9XTXYsYTXY 9301 saLDegITITLEs 200 920094949292989890eLAST ) 
PART 10 COMPUTES 


THE THREE FORCES AND THE THREE MOMENTS AT EACH RING WHICH ARE USED 
FOR THE EQUILIBRIUM CHECKe 


804 


85 


XFOR (K)=20%30e¢14159%ZX( 191) *RA(K)*¥TH(K) 

YFOR(K) =-3014159*XY(291)*RA(K)#TH(K) 

ZFOR(K) =3014159*XY(391)*RA(K)*TH(K) 
XMOM(K)=20*3e1l4159*XYV( 191) ¥*RA(K)*#RACK)*¥TH(K) 

YMOM(K) =3014159*ZX(291)*RA(K) RACK) ¥*TH(K) 

ZMOM(K) =—3014159*ZX(391)¥*RA(K) RACK) *¥TH(K) 

CONTINUE 

PRINT 5559(XFOR(K) sYFOR(K) sZFOR(K) »XMOM(K) »YMOM(K) sZMOM(K) 9 K=198) 


555 FORMAT(///2X7HXFOR(K) SX7HYFOR(K) SX7THZFOR(K 5X 7HXMOM(K)5X 7HYMOM(K) 5 


THIS IS THE END OF PART 10 AND THE START OF PART lle 


1X7HZMOM(K)///(6E12¢5) ) 


END 
PAR ie like THE 


SUBROUTINE GAUSS 2e 
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003920 
003930 
003940 
003950 
003960 
003970 
003980 
003990 
004000 
004010 
004020 
004030 
004040 
004050 
004060 
004070 
004080 
004090 
004100 
004110 
004120 
004130 
004140 
004150 
004160 
004170 
004180 
004190 
004200 
004210 
004220 
004230 
004240 
004250 
004260 
004270 
004280 
004290 
004300 
004310 
004320 
004330 
004340 
004350 
004360 
004370 
004380 
004390 
004400 
004410 
004420 
004430 
004440 
004450 


A OTY Oa 


a 


1 


i 


14 
30 
3.) 
32 
23 


34 
40 


41 


42 
43 


Si) 


THIS IS THE END OF PART 11 AND THE START OF PART 12.6 


SUBROUTINE GAUSS2(NoMsEPoA 9B eX gKER) 
DIMENSION A( 119011) 9Bll91) oX(1191) 
NPM=N4M 


DO 2 KX = loeM 

I = N+K 

DY. 2 J = loeN 
A(JoT) = 8( 59K) 
BDO. 34 b=1ah 
KP=0 

Z2=0-0 


DO 12 K=LoN 

IF (Z-ABSF(A(K0bL)))11912912 
Z=ABSF(A(KoL)) 

KP=K 

CONTINUE 

IF(L-KP)1359209220 

DO 14 J=LyNPM 

Z=A(LoJ) 

A(LsJ)=A(KP oJ) 

A(KP9J)=Z 

TF CABSFCAVL sl) )-EP)50950630 
IF (L-N) 31940940 

LP1l=L+1 

DO 34 K=LP1,N 
IF(A(KoL))32934 932 
RATIO=A(KoeLI/A(LoL) 

DO 33 J=LPl»sNPM 
A(KsJI=A(Ke9J)~RATIO®A(L 9J) 
CONTINUE 

DO 43 I=leN 

IT=Nt+l-] 

DO 43 J=loM 

JPN=J+N 

S=0e0 

IF(TI“N)41943943 

ITIPlL=I +1 

DO 42 K=IIP1l oN 
S=StA( TI eK) *¥X(K oJ) 
XCITTeJV=CAC TI sJPNI-SI/SACIIsI1) 
KER=1] 

RETURN 

KER = 2 

END 


SUBROUTINE DRAW. 


1 
Z 
3 


SUBROUTINE DRAW (NUMPTS» Xe Ys» MODCURV,s 
IL.LLUEs EXSCALEsaYSGALEs 


MODEXAX» MODE YAX » 
LAST) 


A GENERAL CURVE DRAWING AND POINT PLOTTING SUBROUTINE 


PROGRAMMER Je Re WARD 


DATE PEBeml964+s REVISED JUNE 1965 
SYSTEM FORTRAN 60 
OUTe OT LOGICAL TAPE NUMBER 8 


Dh 


PART 2s” Time 


00010 
00020 
09030 
00040 
00c5s0 
00060 
09070 
09080 
99090 
00100 
00110 
00120 
00130 
00140 
00159 
09160 
OO L7G 
00180 
09190 
09200 
00210 
00229 
00230 
00240 
00250 
00260 
00270 
00280 
00290 
00300 
00310 
00320 
00350 
00340 
00350 
00360 
00370 
00380 
00 390 
00400 
00410 
00420 
00430 
004460 
004470 
00000000 
00000010 
00000020 
30 
0 40 
0 30 
0000006 
00000070 
00900080 
00000090 
00000100 


ON VEO OO OOM A OANA NN AD AANA ANN AAA ADOO NAO FOV ON 


ee 


INPUT ARGUMENTS 


le 


2e 


36 


4e 


be 


qe 


Be 


Ge 


NUMP TS 


Piyer 


LABEL 


Digan 


EXSCALE 


YSCALE 


Soret hows. “826 CHS Ne See a es J SCGr 
eles 

ia “NO 3 
eee salts 

eee CF POINTS Seo Seer aie TS MES ALWAYS ae os 


me AT LFAST 25 ANT MUST NOT =XCEEQ) GQmEOume on) iimagamen 
PLOTTING, O9 39 COR CURVE SRAWING. Sages 

cco® 180 
ARRAY OF X-ORQDINATES* DIMENSION AT LEAST EQUAL 02000190 
TO NUMPTS ANN NOT YORE THAN 997 TN CALLING peat Vain 
PROGRAM. 30000210 

S0Gm 220 
ARRAY OF Y-OPDINATES. DIMENSION AS FOR THE Bisa le 
X ARRAY IN THE CALLING PRQOGR24M, 90090240 

000 250 


CONTROES Gige NUNS OF CURVE Ss: ANCOR TS6eTS oF OPE OO ZOG 


POINTSs ON FACH GRAPH, THE CODES AR Caer © 
8, ONLY ONF PLOT ON THIS GRA®OH alae baie itsia 

il FIRS! PLOT ON MULTI —PLOT GRAre wer OO 222) | 

2 INTERMFOTATE PLOT CSN MULTIHWPLOT GFAPHYEAAN3ANN 

2 LAST PLOT ON MULTI=-PLOT GRAPH. OWOGe310 

2060 329 

CONTROLS THE Yep “PUOTs. Iab COUESPARE 00000330 


0 STRAIGHT LINESRIOIN SUCCESSIVE POINTSOGOG 32 eo 


(STANDARD CURVE DRAWING} 90000350 

IL POINTS PLOT LEieWwiit CROSS eX) ODI 360 

Z POINTS PES EDaWh Ta PLUS C+) 00000370 

3 POINTS PLOTTED WITH SSUARE 1000°%384 

4 POINTS PLOTTFD WITH OLTAMOND 90090390 

a POINTS PLOTEEO 21TH TRIANGLE 09009400 

WHEN POINTS ARE BEING PLOTTED {ITYPF=1 THRU. 5)20099N410 
THE POINTS ARE NOT CONNECTED 90000420 
909 430 


IF A CURVE IS BEING DRAWN CITYPF = 0)s LABEL IS 90000440 
A 4=-CHARACTER 8CD CURVE IDENTIFIER WHITH WILL 8&00000450 
REPRODUCED AT THE FND OF THE CURVE. FOR EXAMPLE 300000460 
LABEL = 4H ONE e IF POINTS ARE BEING PLOTTED, 90000470 
LABEL IS AN 8=—CHARACTER [DENTIFIERe THE FIRST & 0900480 
CHARACTERS ARE REPRODUCED WITH THE FIRST PLOTTEDOG 3G sD 


POINTs AND THE LAST 4 WITH THE LAST POINTe SET 90900300 
TO BLANK ANY UNWANTED CHARACTERS. 99000510 

O00 = 520 
AN ARRAY OF TWELVE 8-CHARACTER BCD WORDS. 90009530 
THE FIRST SIX WORDS WILL BE REPRODUCED AS THE ®©)0000540 
FIRS! LINE OF GRAPH TITLE» AND THE LAST SIX 00090550 
WORDS WILL FORM THF SECOND LINEe THE TITLFE MUST 90000560 
INCLUDE THE USERS JOB IDENTIFICATION. DIMENSION 00000570 
12 IN CALLING PROGRAM» AND SET TU BLANK ALL 09000580 
UNWANTED CHARACTFRSe 00000590 

090 600 
X-SCALE (UNITS/INCH) AS POSITIVE FLOATING POINT 00000610 
VARTABLE WITH ONF FIGURE SIGNIFICANCE. SET TO ORONO 15 10) o\ 714, 
ZERO FOR AUTO=SCALE 00000630 

000 640 
Y-SCALE (UNITS/ZINCH) AS POSITIVE FLOATING POINT 00000650 
VARIABLE WITH ONE FIGURF SIGNIFICANCE. SET TO 90090660 


oy) 


ye ae GN GN GG ay) (a Canaan ry Gay OVON YAMA OTN OAA VONAAAANNA 


10. IXuUP 


Vie IYRIGHT 


l12¢ MODE XAX 


132 MODEYAX 


14.6 IWINDE 


156 IHIGH 


166 IGRID 


17. LAST 


NOTE — 


ZERO FOR AUTO-SCALE. 


X-AXIS OFFSET FROM 5SOTTOM OF GRAPH IN INCHES: 
THIS MUST NOT EXCEED THIGHs AND MUST NOT BE 
NEGATIVE. 


SO00C67T0 
$90 680 
000900690 
9C000R a0 
00000710 
COO * M20 


Y-AXIS OFFSET FROM LEFT OF GRAPH IN INCHES. THISO00GO@aG 


MUST NOT EXCFED IWIDEs AND MUST NOT BE NEGATIVE. 


MODE OF X-AXIS OFPFSETs See CODES BELOW. 


MODE OF Y-AXIS OFFSET ss THe GODES ARE ASP FOEEOWS 
0 COMPUTED OFFSET» HOLDING ORIGIN ON 
GRAPHe THE CORRESPONDING IXUP OR 
IYRIGHT IS IGNORED 
l COMPUTED OFFSETs WITH ORIGIN OFF THE 


GRAPH IF APPROPRIATE*s THE CORRESPOND- 


ING IXUP OR IYRIGHT IS IGNOREDe USE 
ONLY WITH AUTO-SCALE 

AXTSSOEFSET GAS cePeGI FLED BY = 1xXUR 70k 
TYRIGHT 


ro 


WIOTH OF GRAPH IN INCHESe THIS MUST NOT EXCEED 
NINEe ZERO WILL BE READ AS EIGHT INCHES. 


HEIGHT OF GRAPH IN INCHESe THIS MUST NOT EXCEED 
SIFTEENe ZERO WILL BE READ AS EIGHT INCHESe 


[IF SET 10 16.4 ONE INCH BY ONE INCH GRID WILE 
BE SUPERIMPOSED ON THE GRAPH. 


INDICATES TO CALLING PROGRAM WHETHER LAST PLOT 
WAS COMPLETED SUCCESSFULLYe THE CODES ARE 


0 EAST PLOT COMPLETED SUCCESSFULLY 
1 PAST PLOT NOT SUGCESSFVE 
2 LAST PLOT NOT SUCCESSFUL AND NO 


FURTHER GRAPH OUTPUT WILL BE ATTEMPT- 


ED UNTIL MODCURV IS NEXT ONE OR ZERO 
5 DRAW WAS ENTERED WITH MODCURV NOT 
EQUAL TO ONE OR ZERO WHILE THE ERROR 
LOCK-OUT WAS SET. 
THIS ARGUMENT MUST ALWAYS BE A NAME IN THE CALL 
STATEMENTe NEVER A NUMBERS 


ALL ARGUMENTS FROM NUMBER 7 THRUs NUMBER 16 ARE IGNORED WHEN 
MODCURV IS EITHER 2 OR 3¢ HOWEVERs ARGUMENTS MUST NEVER BE | 
OMITTED FROM THE CALLING STATEMENTes IT IS MERELY THEIR VALUES 


WHICH ARE THEN 


IRRELEVANT. ARGUMENTS MAY Bt LISTED BY NAME OR 


VALUE IN THE CALL STATEMENTe NO VALUE IN THE CALLING PROGRAM 
WILL BE ALTERED SY THIS SUBROUTINE. 


RFFERENCF --- 


THE BINARY TAPE FORMAT REQUIRED BY THE OFF-LINE PLOTTER IS 
DFSCRIBED IN THE WRITEUP OF THE CDC 160A GRAPH PLOT PROGRAM 
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90900740 
900 750 
9000076N 
000 770 
00909780 
00000790 
00000800 
00000810 
00900820 
00000830 
09009840 
oc0ons5) 
0000086C 
09009870 
000* 880 
00000890 
00000900 
000 910 
00000920 
90000930 
000 940 
00000950 
00000960 
0000 970 
00000980 
00009990 
eXerezab Matera) 
19001010 
00001020 
00001039 
00001040 
00001050 
90001060 
00001070 
60001080 
00001090 
00001100 
00001110 
90001120 
C0001130 
00001149 
00001150 
00001160 
00001170 
00001180 
000011906 
00001200 
90001210 
00001220 


S Cease he Tete FORMAT OFS PVE) Gy THe Cie 6 ee Pew, Sie Se 
a See eT OT MAT OTHE INTE ROCLATION Angee OVS eee ss Vera ty 
& eee SC 
fi eee 2 
IMERIS TONG O)s Y(70Ce -LIITERE(12)e GXTIT(l2)e SYTTI ais eee f. 

1 BET E ECLA )s KAXIS(S1— ICURVIS69O)s GS~IO0C275) 59 TGeN (ees “Geet 28) 

Pa lel Cie ts 2 ) Gna O 
IPOINT = 1 TYPS Pee 39) 

CON ( L203 2 ee ee OC aa 00131) 
Coeeteea= 2 ATT tO Cee ere, [CIR = ClO enya oS Ogee ee. 3 2” 

* Semmes WITH DATS STATEMENT ITM FORTRAN 52-36 O000}1330 
Pee es = | OTN SATA STATEMENT, NONC] 44an 
Poteet = 12506512) QO FIV OC TL ssh 7) VQEO 1] 259 

more’ 1TEST = 0 e260 
JTEST = 73546912 Seas 1° 

99071 CONTINUE CoCr 2 BO 
* REMOVE ABCVE NONSENSE IN FORTRAN 53.4 VO es 70 
G CaEGK PREVIOUS VPERATION OF ROUTINEs (hE VANYseeeCeues ARE 909001400 
C RReST = ¢ IF PREVIOUS GRAPH, IF ANYs COMPLETED 90901410 
(g ITEST = 1 IF PREVICUS GRAPH NCT COMPLETED N0001420 
C ITEST = 2 IF ERROR FOUND WHILE 4OOCURY was CNEs CRaIF “eerie so 
ke MODOC URMeWw A SaatL LEGAL « NOCAL4&49 
Petes = 2b 10%, 10015619C9 09091459 

eo Pre Tt MOOCHRY) 19C3419°%251093 e8001469 
Pome ITEST = C 00901470 
GO TO 1000 00001480 

MOOS TF(MONCURY = 1)10794,1092.21004 90091490 
mo04G LAST = 3 O0001500 
RE TUlan 20001510 

‘@ SET UP ERROR RETURN ROUTINEe ENTRY AT STATEMENT 1005¢ OO001 520 
MOOS IFCITEST)1OC9 +s 1OMEyuCOS 99001530 
1006 IF (MODCURV) 1007219681007 000015490 
moo? PRINT 119C 90901550 
1100 FORMAT ( 59H NO FURTHER GRAPH OUTPUT UNTIL MONCURV NEXT [S ZERO OROD0N91560 
1] ONEse 9/7) 00001570 
Ties) =. 2 00001580 

LAST = 2 90001599 

QE TYRN 90901690 

mo2O8 PRINT 1101 00001619 
11C1 FORMAT (¢ 30H THIS PLOT WILL NOT 3E OUTPUTe 9/) 70001620 
PAST = 1 OVE 1530 
RETURN 00001640 

moOS TFI(MODCURV — 2)101%. 1008 +1010 00001650 
MO10 WetMODCURY — 3):1997s 1Odds1 CC? 000016690 
MO 11 3ITEST = 0 00001670 
60 TO 1098 00001680 

a CHECK LEGALITY @5S INPUT ARGUMENTS e 09001699 
1999 TF(NUMPTS - 2)19292 00001700 
1 PRINT 199 00091710 

199 FORMAT (/5 22H NUMPTS MUST NOT BF LESS THAN 2.6 } 90091720 
GO TO 1095 09001730 

2 LF CIPOINT)9979,9994,90°01 00001740 
9990 PRINT 9170 VOO Cd 20 
SOO FORMAT© (7s 15H ILLEGAL ITYP&e } 09001760 
GO TO 1005 00001770 

SOT PIP CITROINT - 5)9992,9992 «9000 90001780 
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2? 
QNN2 
en) | 


9994 
S005 
9102 


OZ OC 


OZ 1 
HO 2 


1922 
ro2 3 
1924 
ie 3 


Oz 5 


lo 


ey 
18 
ie. 


IF(NUMPTS - 
DRINT 919] 


S03 soa coe 


Mielec; |. S)8 


90991 8N0 

FORMAT (/5 45H “NIVPTS MUST NOT EXCEED 329 FOR OSCTNT PLOTTINGS ) ANAN1E819 
GO TO LACS -0091820 
TF (NEIMOTS = 9NN)2,32,90N6 ANNONA B30 
Ce ieee 2 99001840 
FORMAT (7% 284, NUMPTS MUST NOT EXG2E aro"... ) 90001850 
GO-TO) 1665 00001869 
1X y= 00001870 
IY = lHY 90001880 
AMAXX = -%e2F+109 06001899 
AMAXY = =-0-22F+106 00001900 
AMINX = #922E+1900 NCHOTS To 
AMINY = 49,27£4199 SCAG S20 
90 1€20 T= lsNUMPTS 90001939 
AMAXX = MAX1F(X(1) sAMAXYX) 00001949 
AMAXY = MAXIF(Y( 1) sAMAXY) 00001959 
AMINX = MINIFCXCI) sAMINX) C0001960 
AMINY = MINIFE(YCI1) osAMINY) 00001970 
AMAXA = MAX1LFCABSF(AMAXKX)s ADSFCAMAXY)s ASSFCAMINX)»s ASSFI(AMINY)) 00001980 
IF(AMAXA — 1460£+99)19229102291°21 00001990 
PRINT 1102 00002000 
FORMAT (/»s 58H NO X OR Y VALUF MAY EXCEED 1-0F+99 IN ABSOLUTE MAGNO0092010 
IE OF ABDI A 900002020 
GO TO 1005 00002930 
IF ( ABSF (CAMAXX -— AMINX) ~ 1¢0E&-97)19239102591925 NO0A02940 
IF (ABSF (AMAXY ~ AMINY) -— 1¢NE-97)19249102591025 00002050 
PRINT 1193 00002060 
FORMAT (7/5 38H ALL POINTS HAVF THE SAME COORDINATES. ) 90902070 
GO 10 WeoSs 90002080 
IFCITEST) 49794 00002090 
IF (MODCURV =- 2)5924095 009002100 
IF (MODCURV - 3)6924096 00002110 
PRINT 101 00002120 
FORMAT (/s 17H ILLEGAL MODCURY, ) 900002130 
GO 70,1005 00002140 
I= (MOOCURV) 62998 00002150 
IF CMODCURV = 1169926 C0002160 
FECIWIDE 1051145 1.2 000902179 
Let ¢= 2 Ord 4 Por 00002189 
SPINT Lo2en I Tite Ii 90002190 
FORMAT (/¢9 9H ILLFGAL 9A5929H. GRAPH WILL BE PLOTTED WITH sA5¢ 9090002200 
5H = Be 9/) NOOSZz2Z10 

JWIDE = 8 00002220 
GO TO 14 00002240 
LE GEWIDE. =" 9) 13612520 00002240 
JWIDE = IWIDE 00002250 
IF ( IHIGH) 15516517 00092260 
IT?) = SHIATGH CO092Z270 
PRINT 2G@Zs 1 Tits: tit 00902280 
JHIGH = & 00002290 
GO’ TO 19 90992300 
IF(INYIGH — 15)18518915 90C0Zana 
JHIGH = IHIGH 90002320 
NODFEFXAX = MODEXAX 00002330 
IF (MODE XAX) 20927921 00002340 


wow 


ITIT= B8HMODEKAX, 
PONE JCGsn1TIts Ix 


194 FORMAT (/5 9H [ILLEGAL 9A6&s 32H GRAPH AILL 3 


Also TUAXx = ~~ o/) 
NOOSEXAX = © 
ee TO 27 
IF (MODE KAX -—- 29279722925 
TEC T XP = JU TCH) 249244922 
ITIT = BHIXtHP. 
Sew! LG 6s TEITs [Xx 
NODEXAX = 9 
SO Voe2a 
IF (TXUP8)23426926 
INUP C= 1 XUP 
NODEYAX = MODEYAYX 
Pt 4ODEYAX)28935429 
ITIT=8HMOREYAX. 
PSTN Te 14,6 ITTTs ITY 
NODEYAX = 9 
GOA TO 325 
IF(MODEYAX = £)35430528 
TPC TYRIGHT —- JWI9E)3293203]) 
ITIT = B8HIYRIGHT. 
PRINT 1046 ITITs IY 
NODEYAX = C 
GO TO 35 
Peel YX 1GHT ) 21434434 
JYPIGHT = IYRIGHT 


PT TALIZE PRICR TC SCALING A*!D AXIS LOGATING. 
] POR PASS? FI TH vy RarA. 


IFLAG = 0 FOR PASS WITH XNATAe IFLAG 
DO 2235 IOTA=1,12 
Pot toe whoa = ITITLECIOTA) 


IFLAG = 9 
BFTA = ie 
Senne © EXSCALE 


TAXIS = JYRIGHT 
MODE = NODEYAX 


PSIZE = JWIDE 
Ixy = IX 

LYX = 1Y 

AMAX = AMAXX 
AMIN = AMINX 
SO 10) 5:2 

IFLAG = 1 

BETA = Oe 
SCALE = YSCALE 
TAXIS = JXUP 
MODE = NODEXAX 
ISIZE = JHIGH 
AMAX = AMAXY 
AMIN = AMINY 
™XY = TY 

IY¥X = IX 


CHECK SCALE AND GO TO FIXF OR AUTO SCALE ROUTINES. 


IF ¢SCALFE)53—459456 
PRINT 1145 IXYs IXY 
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ee ne ata 
722902960 
NO 370 
we eS 
Oe. 
Ver AI4AR 
8920241% 
20232420 
900602430 
TON%244 
00602450 
900002466 
20602470 
90002480 
AD902499 
~CoC2 so? 
“OCOZ S15 
OCCeZa20 
OSO92530 
20002549 
COCE2Z 550 
A90902560 
J0GO2Z5'70 
00002580 
09992590 
90002600 
2 O00GZ610 
00002629 
10002630 
90002649 
C0992654 
“0002660 
90002679 
0C002680 
90002690 
00002700 
00002719 
09002720 
90002730 
9002749 
005902750 
00N02 760 
O92 7T7N 
00002780 
00002790 
00002800 
00002810 
90002820 
00002830 
00002840 
00002850 
00002860 
00902870 
090028890 
90002890 
00002900 


ry 


7) 


114 


1 


Biot eyT  C/ » Moi e 
ey lage) 2 = GTS Se 


en 7959 a*agG4 


ryoOvecc Fe ( oe Meee a) Pt Gees (a aurea BO aR EA SS hil fe JSS Wie [= 1 62 See Anna ee 

BE TEISCALE = [QSE4O0NS 7,5 a g8% ~ranznen 
a? re( oe SS bee 1 MF aeSCine Baybee ANAA Dee 
Pee i Ss CA ITS AL = 51 SOA) Soe Sage) AAA 2 
a Re ne ae eee Oe VS ee ee ee A OS ee 
CHECK AND TOMPUTE “8x Ps LOCATION Lo Se ae eee oe 8S ele RAR ADIGA 

CAGE, ITAG © 7 F5 QRIGIS on GRADE Ae 1 Te JT fs sieevcesen, "enema 

fet MOORE al 122 ets te lee ~ AA aes 

193° IITAG =9 2 Oe 
eo The) 2 4 OS ee 
Poenmere iNT Jews 1Y¥Xs TXV¥s I*¥ AAS Se 
1134 FORMAT (/>5 SH YODE241s24H4X 4UST NOT SE 1] UNLESS sAlsS7HSCALEe YS “Ana 
PCATITO=$SCALE )e GRASYH wILE 35 PLEOYTED WIT SUTIC ghls %R=SGA be Cs eee 

nO, aes ' "A °9 ee 

oe Cae CAM AX CS AMIN, = J) .2E—O7 jy IAA a5 O35 eta .999°3— 
eed = SRS AMAY eg HHS 7) LO A451 49,1 NS qg0er2a9gsa 
eae TE (CA YAX VIC 3866917 345 1-377 VOC Se 
Peer AX TS = [Si 2572 “OQC93 is 
Ge To 1C34G J20C3Iie 
PoeemraxlS = 1S1Ze 000024185 
Boy TO 1039 80003140 

Pee TAXIS = 0 OC0C2 tee 
GO-TO 15340 COO0Csie. 

YOS6 IF(SIGNF(lesAMAX! — SIGN (1e94MIN) 119404103951 29460 SLA EOT ORC.) 10) 
1949 ASTZE = ISIZE C0903185 
VAXITS = =AMIN, CAMAX = AMIMYFASTZE +965 “O60 se 

CO 7TO 1030 VA003299 

AUTO SGCALESROUT ENE. 1000 ae 

5° IF(MONE - 1169964959 "00a 
69 AMAX = MAXLE(Nes AMAX) NOOD3 220 
AMIN = MINIF (Qe AMIN) NC9003240 

64 IF(ABSFCAMAX — AMIN) = 1.0€-97)65268568 99093250 
65 PRINT 1169 IX¥9 IXY¥s IYX 00003260 
116 FORMAT (/¢ 5H ALL »Ale47H VALUES EQUAL. AUTO SCALE POSSTGLE ONLY GO0077a 
LIF THE »Al*29H VALUES ARE NON-ZERO AND VODEsAL®s THAX = 2e ) 09903 280 

GO TO) 10s 90003290 
GaeAaSIiZE = TSIZE oOCO4 aa 
SCALE = (AMAX —- AMIN) /SASIZE 0000342 

GO TOe84 OC SNA es 

69 TFE(ARSF(CAMAX - AMIN) - 1¢9£-97) 70974574 09092343 
(OTE (ARSFCAIMAX) = 1. 7F-97)715745,74% NN9QN234N 
fee INT 118s 1TXY 00903350 
118 FORMAT (/s SH ALL »419386H VALUFS ZERO. AUTO SCALE NOT POSSIBLEe )00003360 
GO 719° 1005 190103 32 

74 TF (AMAX = 120E~-97) 76975975 NOCOS3a0 
75 ITFCISIZE = TAXIS) 77s /6977 90003390 
fom>c MLE = 206 90002400 
GO TOs ae 49003410 

fae AK | Se a NNQ032420 
NONE = a 7 = 9009 3a 
SCALE] = AMR TCASIZE = AXIS) QACNA44O 

TR TFECAMIN + J29f —-97) 79979980 NNN See 
TO TE CTAF TS) 8 Tse 9000 Fae 


c 
C 


Bn 
81 
82 
1982 
1382 


} 
2 


1984 
83 
84 
85 
86 


87 
88 


89 
99 


ci! 
92 
ph 
94 
95 


96 
2c 
} 


99 


at 


290 


291 


202 


203 


eget 2See, 


GO, 1D 82 
AXIS = IAX!IS 
SCALE2 = ~AMIN/AXIS 


eG 8 ete 21198465199251994 
PEINT 1IESsMIYXs 1YX 


CON 479 
000090348) 
C9003499 
eT I35C0 
MOOO351C 
SOC C2220 


FORMAT (75 €6H NONE OF THE PLOT LIES Of THE GRAPH WITH THIS) SPEC Treg ra5s5° 


TF sAls&G7H-AXTS LOCATIONS GRAPHY WILL AS PLOTTED WITH YONEs Al, 


THAX =, 9/) 
SODF = ¢ 
SG. | Oy oc 


SCALF = MAXIF(SCELFE1 + SCALE2) 
Bee oC SES TTCSCALEs ISCALIN FACTOR» 3) 
TF(FACTOR -— 5405!385585584 
FACTOR = 1 
ISCALIO = ISCAL1C + } 
GO TO 90 
IF (FACTOR - 2492)87587986 
FACTOR = 5 
GO TO 90 
IF( FACTOR — 14691)8959,5 89588 
FACTOR = ? 
GeaTO 90 
FACTOR = 1 
SeOer = FACTORM@IO.**1SCALIO 
COmMRUTE AXmS LOCATBON IF NECESSARY. AUTO SCALE CASE. 
Be CRORE S= 17925915273 


TF (SIGNF( Le sAM4&X%) = SIGNF (Le sAMIN) 192994992 
TAXIS = -AMIN/SCALE + O45 

ITAG = 0 

Go FO 203 


ITF (AMAX 9959959290 

WAXTS = [SIZE 

SFTA = -AMAX/SCALF 

Ngee ete — Vel+]12)99595596 
PRINT 1299 IXY 


BORMAT (79 15H FHE ORIGIN OF »Als 43H GANNOTeesE OCFFSEI MORE 


ef F+12 INCHES. ) 


GO TO 1005 
IRETA = BETA + 925 
RFTA = =IRFTA 


S574 1S THE NUMBER OF INCHES OF ORIGIN SUPPRESS Ree OsiTive 


TRUE ORIGIN ITS RELOW OF TO LEFT OF THE GRAPH, 
BP GheTA + 16) 97597592 


ITAG = 1 
60 TO 29% 
TAXIS = 9 


BFTA = AMIN/SCALF 

PetoETA = JeE4+12) 26015201996 
HoetTA = FETA + 0.25 

BETA = IBRFTA 

PeCBETA = 16)9332025 202 
PTAG = ] 


RELEASE RESULTS TO REMAINING PART OF PROGRAM. START SECOND 


PASS FOR Y VALUFS IF NOT YFT COMPUTED. 
IF CTFLAG) 20592445295 


61 


SEBS a4 0 
70093550 
DOE 2569 
S0C03i5:70 
00003580 
G2 052790 
92003600 
90003610 
090903620 
206329 30 
90003649 
C090365C 
OLE 266” 
0003679 
60003660 
90003890 
00903700 
00003719 
90003720 
00003730 
00003740 
00003750 
60003760 
00003770 
09003789 
00003790 
90373359 
vigeos 810 
NCGS 829 


THAING 100003830 


99003840 
09003850 
00003860 
90903870 


IF 90003880 


S000389¢ 
NCO9S90C 
00003910 
00903920 
90003930 
00003940 
00903950 
90003960 
00003970 
00003980 
00003990 
90004090 
©0994010 
90094020 


24 


205 


206 


9206 


207 
208 


299 
210 


ra 
22 
aks 


214 
Kaa Ge, 


216 


SEALEY = SCALE 
mX= AGT = FACTSAR 
Poa = 6 [SCALIN 
IXAXIS = IAXIS 
ITAGX = ITAG 
Daizex = ITSIZE 
PETAX = RETA 
GO TO 50 
RETAY = SETA 
SCALEY = SC4LE 
iver.) = FACTOR 
ljvese ln = JSCAL)® 
IYAXIS = TAXIS 
NOW WRITE RECORDS. 
ITAGY = ITAG 
MselzZ£Y = FSi1zZe 
PHP SsGOMPLETESSGASUGULATION Cr SGALE FAGHO {SST Gao 
TAPE (RECORDS. FIRSTs THE SCALE FACTOR TITUES. 
ACT Tt{ 1) =) eet x= 
MATIT(2) = BSHSCALE = 
wert tt 3) = 1[CODStUSGALFEX) 
YATIT(4) = 8H UNITS/I 
epi bitsy = SHNCH, 
Sot tCly = 8H] Y= 
YTI¥(2) = BHSGALE = 
eyTit(3) = ICOMPISCALFY) 
eytit(4) = GH UNITS/1 
JYTIT(5) = 8BHNCH. 
DO 9206 I1=6511 
JIXTIT(I) = 8H 
SY T1T(a) = 38H 


IFCITAGA) 20s 2119207 


IF(BETAX) 20852085209 

wx bIiTt2 ) = 8H ADD - 

GO-TO 210 

JxXtlite) = 2.8 ADD + 

JIxXTT hes) == TeCODr -PSlTAX*®SCALEX) 
JXTIT(9) = 8H UNITS T 
JIXTIT(10)= 8HO ALL X 

Ser itt ll) =SeHVALUES. 


MEU ITAGY) 216% 21642 12 
[EU BETAY) 22352135214 


Ser itty) = 6A ADD - 
GOTO 215 
syEltt?7) = 8H ADD + 
ermine) = [CODE (SETAY*SCALEY) 
eyrirt9) = @H UNITS fF 
ey Tl Pile) =6HO AEE YY 
JYTITCIlI)= BHVALIIFS. 
ICONT(1) = ICONTRL + 4 
INSERT TITLF SIZF (02B) AHEFAN OF MAIN TITLE RECORD. 
SOULS ie oe (ITI PEE s LIIILE) 
TEST FOR ALL BLANK TITLESe 
ICHECK = 8H 
DO 9075 IT=1s6 


9074 3907559N974 
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Pere reel ySTrehean) 


GENERATE 


32994939 
ccoj4nay 


2COgs ee 
/CO9S22o 
{0 Cte 2 
C0904109 
90994119 
‘9C04 123 
00004139 


60094164 
£$9004170 
99009418n 
00004190 
000042909 
C0004219 
90094220 
90004239 
0000904240 
00094250 
00094259 
00094270 
00004289 
00004299 
0000904309 
60004310 
00004 320 
00004330 
00004340 
00004350 
00004360 
09004370 
00004 38N 
90004390 
00004400 
00004410 
00004420 
000044 30 
00004440 
00004450 
00904460 
90004470 
00004480 
90004499 
00004500 
00004519 
000904520 
00004530 
00004540 
00004550 
00004560 
00004570 
00004580 


9974 
9975 


9080 
9081 


9084 
9085 


9090 


309') 


LeCrUPTLe( ly ) 90980s9075e9%80 
CONTINUE 
ITl1 = 1 
PEONT(1) = PE@NT(1]) - 1 
GQ TO 9061 
frit = 0 
DO 99085 I=7512 
IF CITITLECT) = ICHECK) 9984:908529084 
IF CITITLECT))909029085:9099 
CONTINUE 
IT2 = 1 
ICONT(1) = ICONT(1) = 1 
GO 10°9991 
IT2 = 0 
NOW GFNERATE AXES RECORDS, 
LFTMGN = 9*100 
IBOTMGN = 1*100 
Peo= CP TGN 
Jr IBOTMGN +I YAXIS*100 
LH ISIZEX*100 
IHL = LFTMGN + ISIZEX*100 - 107 
KAXIS(1) = IPACK12(1HsJHsLH»s IHL) 
vate = JH — 13 


tt i 


IHL2 = -100 
IVH = (ISIZEX -IXAXIS — L)*IXFACT 
IVH2 = -IXFACT 


KAXIS(2) = IPACKI2(JHLsTHL29IVHsIVH2) 
NH = ISIZEX 

ISH = 8H 14 

IV = LFTMGN + IXAXIS*100 

JV = IBOTMGN 

KAXIS(3) = IPACK1I2(NHs ISHs IVeuV) 

LV = ISIZEY*100 

ivi =" lV = 735 

JVL = IBOTMGN + ISIZEY*100 - 107 
wVEL2 = =-100 


KAXIS(4) = IPACKI2(LVeTIVLsJVL »JVL2) 
TvV = (ISIZEY = IVYAXIS = 1I)*IYFACT 
IVV2 = -IYFACT 

INV = ISIZEY 

ISV = 8H im 


KAXIS(5) = IPACKI2Z(IVVeTIVV2sINVosISV) 
NOW GENERATE CURVES. 

SCX 100¢/SCALEX 

SCY = 100e/SCALEY 

EXAXIS= IXAXIS*#100 

YAXIS = IYAXIS*#100 

ALFTMGN = LFTMGN 


BOTMGN = IBOTMGN 
SHIFTX =EXAXIS - BETAX*®100¢ + ALFTMGN 
SHIFTY = YAXIS - BETAY*100~. + ROTMGN 


ExXSIZE=— 1S) ZEX®100 + LFTMGN + 60 
SIZEX LFTMGN = 60 

NolZe ISIZEY*100 + IBOTMGN + 70 
SIZEY = IBOTMGN = 70 

ICURV(1) = 0 
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49004599 
60904600 
270004610 
201704620 
90904630 
79004640 
09004650 
90004660 
009004670 
00004680 
00004690 
00904700 
90004719 
00004720 
00904739 
00004740 
09004750 
00004760 
00004770 
00004780 
00004790 
00004800 
00004810 
00004820 
00004830 
00004840 
00004850 
000048690 
00004870 
00004880 
00004890 
00004900 
00004910 
00004920 
00004930 
00004940 
00004950 
00004960 
00004970 
00004980 
00004990 
00005000 
00005010 
00005020 
00005030 
00005040 
00005050 
00005060 
00005070 
00005080 
00005090 
00005100 
00005110 
00005120 
00005130 
00005140 


240 
9007 
9700 


9701 
242 


9241 


9242 


241 


9243 


244 
246 
9010 
247 
1260 
260 
261 


Zoo 
9268 


9269 
9270 


9202 


9015 
270 


9020 
Zi2 


273 


IF ( FPOINT )991099007%9010 
IF (XMODF{NUMPTSs2))970099701%9700 
ISWITCH = 1 

GO TO 242 

ISWITCH = 2 

INUM = (NUMPTS + 1)/2 

DO 244 T=l»sTINUM 

Gin= X€2*1-=1)*SCxX. +9 SHIFETX 
Ceea=evite sl —1)*SCY + -SHIETY 
IF (1 -IENUM) 241 99241 9241 

GO TO (92423241) sISWITCH 


3 = Cl] 

C4 = C2 

GO TO 9243 

es X(2*1)*SCX + SHIFTX 


a= Yt2el)*Scy + SHIFTY 
= MINIF(C1sEXSIZE) 
ICl= MAX1F(Cl»s SIZEX) 
C2 = MINIF(C2s YSIZE) 
IC2= MAXIF(C2s SIZEY) 
= MINTFPUGZs9EXSIZE) 
1C3= MAXI1F(C3s SIZEX) 
= MINIF(C4s YSIZE) 
IC4= MAX1F(C4s SIZEY) 
PGURV(CT41) = IPACKI2(91C1ls1C2s1C3s91C4) 
II = INUM + 3 
CALL ITPACKLI(LABEL»s LABEL1s» IDUMMY) 
ICURV(TI-1) = LABEL] 
ICURV(TT) = FCURV4S 
IF ({MODCURV -— 1)2479247%9015 
CALL IREADY (IDUMMY) 
IF ( IDUMMY )500091260295000 
CALL IWRITE (ICONTs IDUMMYs 1) 
IF ( IDUMMY )5000226925600 
CALL IWRITE (JXTITs IDUMMY 911) 
IF ( TDUMMY )50009261 95000 
CALL IWRITE (JYTITs IDUMMYs11) 
IF ( IDUMMY )50009265 5000 
Her Cl1T1)926949265859269 
CALL IWRITE (LTITLEs IDUMMY»s 7) 
IF ( IDUMMY )500099269 95000 
VEC 2 S27 159270592 71 
BALt IWRITE CLTATleEt(7)s IDUMMY, 7) 
IF ( IDUMMY )500099271s5000 
CALL IWRITE (KAXISs IDUMMY>» 5) 
IF ¢« IDUMMY )5000s9015%5000 
IF (IPOINT)9020927029020 
CALL IWRITE (ICURVs IDUMMY>s IT) 
IF ( IDUMMY ) 5000990205000 
HE (MODCURV — 1)272527239025 
me (IGRID — 9 1)9025527399025 
GENERATE GRID IF CALLED FORe 


IX100 = ISIZEX*190 
Iy¥TOO ="91S1ZEY#in0 
NEXT1 = IBOTMGN 

NEXT2 = LFTMGN + IX100 
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00005150 
00005160 
90005170 
00005180 
00005190 
09005200 
00095210 
00005220 
00005230 
00005240 
00005250 
00005260 
00005270 
00005280 
00005290 
00005300 
00005310 
00005320 
00005330 
00005340 
00005350 
00005360 
00005370 
00005380 
00005390 
00005400 
00005410 
00005420 
00005430 
00005440 
00005450 
00005460 
00005470 
00005480 
00005490 
00005500 
00005510 
00005520 
00005530 
00005540 
00005550 
00005560 
00005570 
00005580 
00005590 
00005600 
00005610 
00005620 
00005630 
00005640 
00005650 
00005660 
00005670 
00005680 
00005690 
00005700 


sg al ie 


1274 
ie 15 
1276 


1277 


1278 


1279 
1289 
1281 


90:25 


9030 


9031 
9032 
9033 
9034 


9635 


9036 


9037 


9038 


JGRID@) = 0 
DO 1274 J=191992 
JGRIDO(J+1) = IPACK12 (LFTMGNs NEXT1ls NEXT2s NEXT1L) 
IF(NEXT1 - IBOTMGN - 1Y100)12729127551275 
NEXT1 = NEXT1 + 100 
JGRID(J+2) = IPACK12 (NEXT2s NEXT1l» LFTMGNs NEXT1) 
IF(NEXT1 = IBOTYGN — 1Y100)1274%127691276 
NEXT1 = NEXT] + 100 
JGURID( J+2) IPACK12 (NEXT2s NEXT1s NEXT2s NEXT1) 
JGRID(J+3) ICURV3 
JGURID(J+4) ICURVS 
CALL IWRITE(JGRID»s IDUMMYs J+4) 
IF ( IDUMMY )590091277+5900 
NEXT1 = LFTMGN 
NEXT2 = IBOTMGN + IY100 
DO 1279 J=lells2 
JGRID(J+1) = IPACK12 (NEXT1s IBOTMGNs NEXTls NEXT2) 
IF(NEXT1 — LFTMGN = 1X100)1278,128051280 
NEXT1 = NEXT1 + 100 
JGRIO(J+2) = IPACK12 (NEXT1»s NEXT2s NEXT1l»s IBOTMGN) 
IF(NEXT1 — LFTMGN = 1X100)12794,128191281 
NEXT1 = NEXT1 + 190 
JGRID(J+2) = IPACK12 (NEXT1s NFEXT29 NEXTls NEXT2) 
JGRID( J+3) ICURV3 
JGRID( J+4) ICURV4S 
CALL IWRITE (JGRID»s IDUMMYs J+4) 
IF ( IDUMMY )50002*9025%5000 
IF ( IPOINT )9030927629030 
GENERATE POINT PLOT RECORDS IF CALLED FOR. 
1OuT =_0 
CALL IPACKL1 (LABEL» LABEL1]:»> LABEL2) 
DO 9050 T=1sNUMPTS 
Cl = XC 1)*SCX + SHIFTX 
C2 = Y(I)*SCY + SHIFTY 
IF(Cl — EXSIZE)903199031%9034 
IF(C2 = YSIZE)903229903239034 
IF(Cl =— SIZEX)9034 9903399033 
IF(C2 — SIZEY)90349903529035 
IOUT = IOUT +1 
GO TO 9050 
Icl = ¢] 
GO TO (90363s9037s9038%9039 99040) sIPOINT 
GENERATE CROSSe 
ICURV(2) = IPACK12 (1C1-59 I1C2-5s9 I1C1+55 I1C2+5) 
ICURV(3) = IPACK12 (ICl » IC2 9» I1C1-5e9 +€2+5) 
ICURV(4) = IPACK12 (101459 102-59 I10C1+55 IC2-5) 
GO TO 9041 
GENERATE PLUSe 
ICURV(2) = IPACK12 (IC1l +» I1C2=5e I1C1l » IC2+5) 
ICURV(3) = IPACK12 (IC1l 2» IC2 » IC1-5s I1C2  ) 
ICURV(4) = IPACK12 (1C145%s I1C2 9» I1C1+5s IC2 ) 
GO TO 9041 
GENERATE SQUARE. 


ICURV(2) = IPACK12 (1C1+4s I1C2-45 I1C1+4,s I1C2+4) 
ICURV(3) = IPACK12 (1C1-4s 102449 IC1-4s IC2-4) 
ICURV(4) = IPACK12 (IC1+4s I1C2-459 IC1+4s I1C2-4) 
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00005710 
00005720 
00005730 
00005740 
00005750 
00005760 
00005770 
00005780 
00005790 
00005800 
00005810 
00005820 
90005830 
090005840 
00005850 
00005860 
00005870 
00005880 
00005890 
00005900 
00005910 
00005920 
00005930 
00005940 
00005950 
00005960 
00005970 
00005980 
00005990 
00006000 
00006010 
00006020 
00006030 
00006040 
00006050 
00006060 
00006070 
00006080 
00006090 
00006100 
00006110 
00006130 
00006140 
00006150 
00006160 
00006170 
00006180 
00006190 
00006200 
00006210 
00006220 
00006230 
00006240 
00006250 
00006260 
00006270 


* oY 


9939 


9040 
904] 
9042 


9043 
9044 


9045 
9046 
9050 


99048 
9104 


276 
PAA | 
278 


Ie 1% 


656 


279 


5000 
5001 


5002 
5100 


5670 


GO TO 9041 
GENERATE DIAMOND. 


ICURV(2) = IPACKI2 (1C1+452 1G? Weel Cl. gules 5) 
PCURV(3) = IPACKT2™°(TC1=Ss™ Games 161 wamre2e—5) 
TCURV(4) = IPACKI12 (IC14#5s [C2 9 1¢C14+55 [c2 
GO TO 9041 

GENERATE TRIANGLEe 
ICURV(2) = IPACKi2 (1C1+5s IC2=-3s IC1l » 1062+6) 
ICURV(3) = IPACK12 (1C1-5s I[C2-35 I1C1+5s I1C2-3) 
ICURV(4) = ICURV(3) 
IF (IT — NUMPTS)9043 9994239043 
ICURV(3) = LABEL2 
GO TO 9046 
IF(I = 1)90459904459945 
ICURV(5) = LABEL1 
GO TO 9046 
ICURV(5) = ICURV3 
ICURV(6) = ICURV4 
CALL IWRITE (ICURVs IDUMMY>s 6) 


IF ( IDUMMY 15000 995055000 
CONTINUE 
IF (10UT)99048927639048 
PRINT 91049 JOUT 
FORMAT (/s 1Xo I2s 
SET UP RETURN. 
IF (MODCURV) 27792789277 
IF(MODCURV - 3)27992789279 
PLESt 2=-70 
PRINT 130e(JJTITLE(I)sTI=1l912) 
FORMAT (/s 19H GRAPH TITLED e . 96A8s/919Xs6AB8» 
24H ¢« e HAS BEEN PLOTTEDe 29/91H0) 
IDUMMY = ITYP2¢1DUMMY) 
IF ( I DUMMY )56703656s5670 
LAST, = 0 
RETURN 
ITEST = 1 
IDUMMY = 
LAST = 0 
RETURN 
THESE ARE THE NORMAL RETURNSe 
NOW SET UP THE RETURN FOLLOWING A TAPE ERROR. 
IF (MODCURV - 1)5001%500125002 


29H POINT(S) WERE OFF THE GRAPHe 39/) 


ICLOCK (I DUMMY ) 


IDUMMY = ITYPE1(IDUMMY) 

GO TO 247 

PRINT 5100 

FORMAT (/s 36H TAPE ERROR IN WRITING GRAPH OUTPUTe ) 
IDUMMY = ITYPE1(IDUMMY) 

GO TO 1607 

IDUMMY = ITYPE1(IDUMMY ) 

END 


SUBROUTINE IREADY fT timmy) 
SELECTS TAPE 8 (. LL LOOP UNTIL READY)». WRITES EOF ON 8e 
MACHINE LANGUAGE wILL NOT BE NECESSARY IN FORTRAN 62-36 
LOC(IFIVE = 5)e 
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00006280 
00906290 
09006300 
09006310 
00006320 
00006330 
00006340 
00006350 
00006360 
00006370 
00006380 
00006390 
00006400 
00006410 
00006420 
00006430 
00006440 
00006450 
00006460 
00006470 
00006480 
00006490 
00006500 
00006510 
00006520 
00006530 
00006540 
00006550 
00006560 
00006570 
00006580 
00006590 
00006600 
00006610 
00006620 
00006630 
00006640 
00006650 
00006660 
00006670 
00006680 
00006690 
00006700 
00006710 
00006720 
00006730 
00006740 
00006750 
00006760 
00006770 
00006780 
00006790 
00006800 
00006810 
00006820 
00006830 


FXF (529413) Exe O20 0 5:) » SELECT READ AND WAIT TAPE. AC006840 

YWEX FXF7T(0995%R) Soe LADY je iY 7 CONS CA eACT | Vitae OO90S5859 
PAR IU 529° 7B) SLU CTNEX dee EXIT ON TAPE READY 0006869 

BOA (IFIVE) SAU CTSUFO. TFRAINATE 09906870 

18U= EXF5(N)« BUFFER. 00006880 
YRDY EXF (520418) EXF7(520008). SELECT AND WAIT TAPE 8-e 00006890 
ENA (O)e CLEAR Ae 00096909 

EXF (020008) EXF (520063). STGP CLOCK AND BACKSPACE 8e 00006910 

EXF7( 520013) SLI CIEND). EXIT IF NOT AT LOAD POINT. 00006920 
~EXF7(520003)6 WAIT TAPE 8e 000069 30 
EXF7(520078) SLI (1LEOF). EXIT IFeNO GOR, 00006940 

ENA (IDUIMMY) SAU (2BUF).e MOVE 00006950 

INA (1) SiS Gey E ls FORWARD 00006960 

2BUF EXF5(N) EXF7(5200CB). OVER RECORD. 00006970 
TEOF ENA (0) EXF7(00C618B)- CLEAR As WAIT CH 66 00006980 
EXF (620418) EXF7 (620093). SELECT WRITE AND WAIT TAPE. 00006990 

FXF (620418) EXF7(62090B). SELECT AND WAIT TAPE 8e 09097900 

EXF (620038) EXF7(620008). WRITE EOF AND WAIT. 00007010 

EXF7 (620078) ENA(10)~. EXIT ON NO END OF TAPE. 00007020 

lEND STA CIDUMMY)e 00007030 
END 00007040 

(a 00007050 
SUBROUTINE IWRITECISTART» IDUMMYs IWORDS) 00007060 

e WRITE RECORD OF IWORDS» STARTING WITH ISTARTse PUT IDUMMY = 00007070 
C TEP RECORD CORRECTLY WRITTENs OTHERWISE SET NON-ZERO.’ 00007080 
% MACHINE LANGUAGE WILL NOT RE NECESSARY IN FORTRAN 62-36 00007090 
LOCC ISIX = 6)e 00007100 

~EAie fOC0618)) 6 WAIT CH 66 00007110 

EXF (62041B) EXF7(62000B). SELECT WRITE» WAIT TAPE. 00007120 

FXF (62041fA) EXF7(62090B). SFLECT AND WAIT TAPE 8e 00007130 

ENG (111B). SET COUNTER. 00007140 

1AGN ENA (ISTART) INAC1)¢ STARTING ADDRESS» 00007150 
SAL (1BUF) ADD( IWORDS). TERMINAL ADDRESS» 00007160 

1BUF STACISIX) EXF6(N)« BUFFER OUTe 00007170 
ENA(0O) EXF7(62000B). CLEAR Ae WAIT TAPE 8e 00007180 
EXF7(62007B) SLU (TEND). EXIT IF NO END OF TAPEe 00007190 
FXF7(62093B) SLI (2AGN)6 EXIT IF NO PARITY FRRORe 00007200 
EXF7(620948) SLU (2END).« EXIT IF LENGTH ERROR. 00007210 

2AGN EFXF (62006R) EXF7(620090B). BACKSPACE AND WAIT. 00007220 
QRS (3) QJP1(1AGN). TRY WRITE 3 TIMESe 00097230 

emt EXF (620038) ENA (10). WRITF EOFs NON ZERO Ae 00007240 
ZEND STA (IDUMMY )e STORE RESPONSEe 00007250 
END 00007260 
90007270 

FUNCTION ITYP2 CIDUMMY ) 00007280 

TYPE WORD GRAPH. 00007290 

WILL NEED REWRITING IN FORTRAN 62-36 00007300 

CON(LC = 57Be M1 = 45131230150N50C00B )e 00007310 
LOC(ITWO = 2)e 00007320 
~EXF7(O00061B). WAIT CH 66 00007330 

FXF (620419) EXE 70620003). SELECT AND WAIT TAPF. 00007340 

FXF (629418) EXF7(62000B), SELECT AND WAIT TAPE 8e 90907350 

EXP MWe 201038) EXF (0199008). WRITE EOFe START CLOCK. 00007360 

ENA (0) EXF7(62000B). CLEAR Ae WAIT TAPE 8e 00007370 
FEXF7(620078) ENA (10) EXIT IF NO END OF TAPEes 00007380 

STA (ITYP2)e STORE RESPONSE. 00007390 
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KAN 


1TYy°? 


Die 
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[EXE (COOOZTBR) 6 WAIT CHigz. 
EXF7(11141B) SE eT YP i. EXIT IF UPPER CASE. 
EXF 4211006) ENA (LC+l)e TYEE 
STA (ITWO) EXF2(LC)e LOWER CASE.e 

-EXF7(000218). WAIT CH 2e 
EXF (211008) ENA (Ml+l)e TYPE 
STA (ITWO) EXF2(M1l)« GRAPH 

EN) 

FUNCTION ITYPE1] (1DUMMY ) 


REWIND TAPE 8s REQUEST NEW TAPEs AND WAIT TILL READY. 
WILL NEED REWRITING IN FORTRAN 62-36 


CON(LC = 57Bs Ml = 4315112030242004B»s M2 = 1103302204062031Bs 
M3 = 0401301520043342B). 

RSV(MESS = 3). 

LOC(ITWO = 2)e 
-EXF7(00061B). WAIT CH 66 

EXF (620418) EXF7(62000B). SELECT AND WAIT TAPEe 

EXF (620418) EXF7(62000B). SELECT AND WAIT TAPE 8e 

EXF (62003B) EXF 7 veZ0098)', WRITE EOF AND WAIT. 

EXF {(62007B) EXF7{(00021B). REWIND WITH INTERLOCKs WAIT CH 
EXF7(11141B) Sia ClTY Pi EXIT? [FP °UPPERBGASE. 

EXF (211906) ENA (LC+1l)e TYPE 

STA (ITWO) EXF 2ULG)% LOWER CASEe 

LDA (M1) STA (MESS). 

LDA (M2) SIA (MESS+I). 

LDA (M3) STA (MESS#2). 
=EXR? (O0021B jes WAIT CH 2e 

EXF (21100B) ENA (MESS+3). TYPE 

STA (ITWO) EXFZ(MESS).« MESSAGE. 
~EXF7(00061B) . WAIT CH 66 
—-EXF7(62000B) WAIT TAPE. 


EXF (62041B) EXF7(62000B). SELECT AND WAIT TAPE 8e 
EXF (010008). START CLOCKe 
END 


FUNCTION ICODE (ANUMBER) 
CODES ABSOLUTE VALUE OF A FLOATING POINT NUMBER (BETWEEN 
1¢0E-100 AND 1-0E+100) INTO 8=-CHARACTER BCD WORD OF THE FORM 
1¢623E+45~¢ ICODE = 8HO.O0E+00 IF MAGNITUDE OUT OF RANGE. 
CHECK AVAILABILITY OF LIBRARY FUNCTIONS IN FORTRAN 62-36 

DIMENSION II(8) 

BNUMBER = ABSF(ANUMBER) 

IF(BNUMBER — 120E+100) 79696 

IF(BNUMBER - 1.20E£-100) 69692 


ICODE = 8HO.00E+00 
RETURN 
WHIS 1S ERROR@EXITs 


CALL SCALEIT (BNUMBER» ISCAL10»5 FACTOR» 3) 


ISIGNSC = XSIGNF(1sISCALI10) 
ISCAL10 = XABSF(ISCAL10) 
IFACT = FACTOR*100-001 


I1(8) = XMODFC(ISCAL10910) 
II {(7)=I1SCAL10/10 

IF( ISIGNSC) 49393 

L1(6) = 8H + 
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00007400 
00007410 
00007420 
00007430 
00007440 
00007450 
00007460 
00007479 
00007480 
90007490 
00007500 
00007510 
00007520 
00007530 
00007540 
90007550 
20007560 
00007570 
00007580 
00007590 


200007600 


00007610 
00007620 
00007630 
00007640 
00007650 
00007660 
00007670 
00007680 
00007690 
00007700 
00007710 
00007720 
00007730 
00007740 
00007750 
00007760 
00007770 
00007780 
00007790 
00007800 
00007810 
00007820 
00007830 
00007840 
00007850 
00007860 
00007870 
00007880 
00007890 
00007900 
00007910 
00007920 
00007930 
00007940 
00007950 


ay 


ROAOADN 


xa a) 


An 


10 
11 


12 


13 


INEX 


€ 2) 

Lict6)* =.8H = 

Prt)? =a0r 2 

1104) = XMODF(TIFACTs:9) 

PIt2) = (CAROCMTTPATT s 100) 770 

11(2) = 8H e 

Pi(1) = )EASi7 19s 

CALL ITPACK (IIs IPACKEHD) 

ICODE = TPACKEN 

RFTURN 

FND 

Sean OUT) 
FINDS FACTOR (BETWEEN 1eDd AND Ge9Fe ee) 
DEE! 


BE 6B 


ISCALLO=LOGLOF (ANUMBER) 
FACTOR = ANUMBER/1Ce*## I SCALIG 
PEtCEACTOR -— Cel)le2s2 

FACTCR = FACTOR*1900. 

ISCAL10O = ISCAL19 = 2 

GO TO 8 

IFCFACTOR - 140)398s4 

FACTOR = FACTOR*10¢ 

foeAelO = ISCALI10O - 1 

GO TO 8 

IF (FACTOR = 100¢9)69595 
FaCTOR = FACTOR/1006 

ISCAL10 = ISCAL10 + 2 

GO TO 8 

IF (FACTOR - 10¢0)89797 

FACTOR = FACTOR/106¢ 

ISCAL1O = ISCAL10 + 1 

IF (MODE )9%9910 

MODE = 6 

ec TO 11 

IF (MODE — 10)1191199 

IFACTOR = FACTOR*10¢e** (MODE - 
FACTOR = IFACTOR 

FACTOR = FACTOR/10e**( MODE - 
IF (FACTOR -— 10¢)13912912 
FACTOR = le 

ISCALIC = ISCAL10 + 1 

RETURN 

END 

SUBROUTINE IPACK (IIs IPACKED) 


TAKES 8 SIX-BIT WORDS AND PACKS THEM LEFT TO RIGHT 
PACKEDe IF WORD IS ZERO» 
CONVERT TO CODAP FOR FORTRAN 


TT 


NE SCALEIT (CANUMSER » ISCAL19 s FACTOR » MODE) 

AND SCALE OF 10 AS 
NED BY ANUMBER = FACTOR*®1Ce®*ISTCALIO® 

MODE IS™THE NUMBER Or Sieh Fr ICANT 
will ape PUT 
CHECK AMATWABILITY.OF LOGROF. IN FORTIAN G2=36 


ETWEEN 1 AND 10 CR IT 


CON(TIZERO = 12B)-6 


SHUI PSAVvE) 
LDAL( IT) 
LDA (IZERO). 


ENI1(8)-6 
AJP1(2NEX)eo 


+ 025 


12B 1S SUBSTITUTED. 
62-36 
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SQuAt 19 SIXs 


ICOC7TSBO 
~00979350 
70908900 
00008010 
COO %ea20 
00008930 
%0908940 
“9008050 
00008060 
00008070 
00008080 
09008090 
90008109 
%C008119 
20006120 
00008130 
00008140 
90008150 
60008160 
90908170 
09008180 
00008190 
90008200 
00008210 
00008220 
00008230 
C0008240 
00008250 
00008260 
000608270 
00008289 
90008299 
00008300 
00008310 
00008320 
00008330 
00008340 
C0008350 
00008360 
00008370 
00008380 
00008390 
C0908400 
00908410 
00098420 
00008430 
00008440 
00008450 
00008460 
00008470 
00008480 
00008490 
09008500 
00098510 


2NEX 


* OY 


1NEX 


* 


[Roe os) 
ESkit= 1) 


STO ( TPACKED) 


END 


SUBROUTINE ISHIFT6 


INT 1(=2) 6 


SLU (1INEX)e 
ISLULCISAVE) « 


C171 DWE» 


Little 


INSERTS 02B AHEAD OF 6-WORD TITLE RECIRDe 
WILL HAVE TO BE CONVERTED TO CODAP IN FORTRAN 62-36 


WATCH ARRAY 


CONC IBLANK = 
STU1(ISAVE) 
ENA (2)e 
EDG@I(CITITLE) 
STALCLIITEE) 
ISK1(6) 

LDQ (IBLANK) 
ENT1(7) 
ENA (2)e 


EDO RCW IWEE ) 


PNT 1g) 
Nd C=) 
TSKIC12) 


LDQ (IBLANK) 


ENTI1(14) 


LIUL(CISAVE de 


END 


FUNCTION IPACK12 


ENT 1(1)e 


LLS (42).6 
LLS (6)e 


SLJ CINEX)- 


LLS (42)6 


SHATCLIITLE)« 


LLS (42)e¢ 


>) ASCE EE)» 


LLS (6)e 


SLU (2NEX)e 


LLS (42)e6 


STAVTOETIILE) « 


(TONE sI2sI139T4) 


INDEXING IN FORTRAN 62-36 
2929202020202C029R) 


SAVE INDEXs SET COUNT. 
ENTER O2Be 
PERFORM 
SHIFTINGe 
CHECK TReCOMPEETEs 
COMPLETE LAST WORD. 
STORE LAST. 
REPEAT 
FOR 
SECOND 
WeTEE 
LINEe 


RESTORE INDEX. 


PACKS FOUR 12-BIT WORDS INTO ONE 48~-BIT WORD. 
WILL REQUIRE CONVERSION TO CODAP IN FORTRAN 62-36 


LDA( TONE) 


QLS(36) 


LDQ( 13) 
Est 12) 


QOLS(36) 
STAC TPACK12)e 


END 


SUBROUTINE TPACKLI1 
PACKS TWO 4-CHARACTER LABFLSe 


LDQ( 12). 
LUS(l2ie 
QLS(36). 
LDQ( 14). 
BES C12). 


(LABELs 


LARFL1s LABEL2) 


USE DECODE/ENCODE IN FORTRAN 62-36 
CONC IFLAG = 37773777B)e 


LLS(24) 
LDA (IFLAG) 
LDA (IFLAG) 


STA (LABEL2)-« 


END 


STA (LABEL1)-. 
LDQ (LABEL). 


LLS(24)e 


FUNCTION ICLOCK(IDUMMY ) 


FXF (O1000B). 


END 
lee = C2 
END 
END 


START CLOCK. 
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€9008529 
C09008530 
00008540 
60098550 
90008560 
00008570 
00008580 
00008590 
00008600 
00008610 
09008620 
00008630 
00008640 
00008650 
00008660 
00008670 
00008680 
00008690 
00008700 
00008710 
00008720 
00008730 
00008740 
00008750 
00008760 
00008770 
00008780 
00008790 
00008800 
00008810 
00008820 
00008830 
00008840 
00008850 
00008860 
00008870 
00008880 
00008890 
00008900 
00008910 
00008920 
00008930 
00008950 
00008940 
00008960 
00008970 
00008980 
00008990 
00009000 
00009010 
00009020 
00006120 
00009030 


APPENDIX B 
PROGRAM NOMENCLATURE 
The following listing of program symbols is provided to help the 
reader understand symbols used in the digital computer program, Symbols 
identified in section 3 as well as symbols readily identifiable are not 


listed in this appendix, 


X(Y) individual strain gage readings 

¥(I) percentage of depth of individual strain gage readings 
AZ the zero depth intercept of strain axis 

ANG(L) the angle of Ome from gage leg number 1 

SMAX(L) Jima x | 

SMIN(L) Tmin 


THETA(L)  ©p 
SIGXX(LL) experimental value of Gy, 


TRUXY(L) experimental value of Txy 


PHI(L) the angle ¢g 
XSXX(L) components of experimental values of G, in proper form 
YSXX(L) for plotting 
XTXY(L) i components of experimental values of xy in proper form 
XTXY(L) for plotting 
ZX(1,2) 

1 seven coefficients used in computing smoothed value of On 
ZX(7,1) 
X¥(1,1) | 

seven coefficients used in computing smoothed value of Txy 

XY(7,1) 
ZZZ the angle in radians on the circumference of the piping between 


two consecutive smoothed values of Qn or Txy 


PHI 1(1) the angle in radians from rosette number 12 of a smoothed 
value of Gy orTky 


SIGN 7(I) smoothed value of Gn. 
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TAUN 7(I) smoothed value of Txy 


XSNR(T) ; ; 
YSNR(I) components for plotting zero stress circle on On plot 
XTNR(T) 


YINR(I) components for plotting zero stress circle on Txy plot 


YSN7(I) plotting 


XTN7(I) 


components of smoothed value of Txy in proper form for 
YTN7(T) 


XSN7(I) i components of smoothed value of Gnu in proper form for 
| plotting 
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APPENDIX C 
DESCRIPTION OF SUBROUTINE GAUSS 2 


A. IDENTIFICATION 
“TITLE: Solution of Simultaneous Linear Algebraic Equations 


CO-OP ID: F2 

FORTRAN: Gauss 2 

CATEGORY: Simultaneous Linear Equations 
PROGRAMMER: CC. B. Bailey and Mary Haynes 
DATE OF COMPLETION: April 1963 


B. PURPOSE: To solve one or more sets of linear algebraic equations using 
Gaussian elimination with row pivoting and back substitution. 


G. USAGE: 
1, Calling sequence 
Call Gauss 2 (N, M, EP, A, B, X, KER) 
N - order of the matrix in the equation Ax = b, that is, the 
number of linear equations (maximum = 50). 
M = number of vectors b for which solution vectors x are to be 
obtained, that is, the mumber of sets of linear equations (Maximum = 60), 
EP ~ matrix condition parameter (see Mathematical Method). 
A - the elements of the matrix of coefficients of the equations, 
stored in the form (Ap 
t= 1, 2, cosy N 
BT po? Ss uveiajers, 
A is dimensional as a 50 x 50 array. 
B = the elements of the column vectors b etc., stored in the 
form (B, 5) 
' 2=1, 2, ..., 8 
J=1, 2, ooo, M 
B is dimensional as a 50 x 60 array. 


X - the components of the x= vectors ‘stored in the form (X44). 
2 a aie a. 
j=1, 2, ..,M 
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X is dimensional as a 50 x 60 array, 


D. 


36 
he 


De 


KER - error flag (See Error Return) 
Space Required: 


Temporary or Common Storage Required 


Error Return: Argument KER 


1. indicates no errors 
2. indicates that matrix is singular or nearly singular, 


Accuracy: Not applicable. 


METHOD 


This subroutine is an adaptation of the program F2 VTEX LINEQN, 


Please consult the write-up of that program for the mathematical method. 


The subroutine gives no printout. 


Th 


APPENDIX D 
DESCRIPTION OF SUBROUTINE DRAW 

A. IDENTIFICATION: 

TITLE: General Graph Output Subroutine 

CO-OP ID: J7-NPS~DRAW 

CATSGORY: Output for Off-Line Plotting 

PROGRAMMER: J. R, Ward 

DATE: February 1964; REVISED JUNE 1965 
B. PURPOSE: 

This subroutine, when provided with the necessary information, generates 
a magnetic tape in the proper format for subsequent off-line graph plotting 
using the CDC 160 or 160A GRAFPLOT program and a CalComp 165 Plotter 
(see references 1 and 2), Provision is made for curve drawing and point 
plotting, automatic scaling, graph titling and axis annotating. An 
attempt was made to provide a considerable amount of flexibility, at the 
expense, necessarily, of a relatively large number of arguments and a 
rather high memory requirement. 
CG. USAGE: 

1. Definitions: 

In what follows the word "graph" will be taken to mean one piece 
or frame of graph paper on which there may be plotted one or more curves 
and/or sets of points. A “curve” will mean a, contimuous line generated 
by the sequence of straight lines joining successive points of the best 
defining the curve. A “point plot™ will describe the representation of 
a succession of points by means of symbols (such as a cross) on the graph. 
The points are not connected in a point plot, 

2. Calling Arguments: 


All necessary information is transferred to DRAW through the calling 
15 


arguments. The call statement is: CALL DRAW (NUMPTS,X,Y,MODCURV, 


ITYPE, LABEL, ITITLE, EXSCALE, YSCALE, IXUP, IYRIGHT ,MODEXAX , MODEYAX, IWIDE, 


IHIGH, IGRID, LAST) 


It is important to realize that one and only one curve or set of points 


is plotted each time DRAW is called. However, it is possible to call 


DRAW repeatedly if several curves and/or sets of points are wanted on 


one graph. 


The calling arguments are as follows: 


ae 


De 


d. 


NUMPTS : 


MODCURVs 


The number of points defining a curve (2 S NUMPTS 

< 900), or the number of points to be point plotted 
(2 < NUMPTS < 30). 

The array of X-ordinates ( ]X,| £1079 for L=ly 2 seen 
NUMPTS). X must be dimensioned at least equal to 
NUMPTS in the calling program. 

All points will be considered to have the same 
X-ordinate if |X... - Xnin| < 10777. The common 
value will be put equal to zero if eater | S$ 10-97, 
The array of Y-ordinates, with properties corresponding 
to the.X-ordinates, above. Y must be dimensioned at 
least equal to NUMPTS in the calling progran. 

Controls the number of curves, and/or sets of points 

on one graph: 

This is the only curve, or set of points, to be plotted 
on this graph. 

This is the first of two or more curves, and/or sets of 
points, to be plotted on this graph. 

This is an intermediate curve, or set of points. 


This is the last curve, or set of points, for this graph. 


76 


Co 


f. 


ITYPE: 


LABEL: 


tl 
Wo 


=4 
=D 


Controls the type of plot (i.e., curve or point plot): 
This set of points is to be represented by 4 curve, 
These points are to be plotted with a cross ee. 
These points are to be plotted with a plus (+). 
These points are to be plotted with a square (C0 ). 
These points are to be plotted with a diamond (O ), 
These points are to be plotted with a triangle ( A ), 
This is a Hollerith curve or point identifier. Ifa 
curve is being drawn, LABEL mst have 4 characters 
(including any blanks), and these will be reproduced 
beside the end of the curve. This argument can be set 
in the calling program by a statement such as 

LABEL = 4H ONE (= blank) 
or LABEL = 4H1234 


or LABEL = 4H The latter must be used when 


AAAA 
no label is wanted. 
If a set of points is being plotted, LABEL is an 8&- 
character identifier, The first 4 characters will be 
reproduced beside the first point, and the last four 
characters will be reproduced alongside the last point. 
This argument can be set by statements such as 

LABEL = 8HFRSTLAST 
or LABEL = 8H,,,, ONE 


or LABEL = 8H,ONE,123 


or LABEL = 8H The latter must be used if 


AAAAAAAA 
no labels are wanted, 
The above arguments, a. through f. (and q.), have 


meaning every time DRAW is called. On the other hand, 
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Be 


ITITLE: 


EXSCALE: 


YSCALE: 


IXUP: 


IYRIGHT: 


MODEXAX:s 


the remaining arguments, g. through p., have no 
meaning except when MODCURV = 0 orl. 

An array of twelve 8-character Hollerith words, the 
first six of which will form the first title line, and 
the last six the second. The array must be dimensioned 


12 in the calling program, must contain the user's job 





identification, and must have unwanted characters set 
to blank. For example: 

DO1lI=1,12 
1 ITITLE(I) = 8H 

ITITLE(1) = GH, SMITH, ,, 

ITITLE(2) = 8HJe Jenny 

ITITLE(7) = 8H,TESTIT, 
X-scale in units per inch (10777 < EXSCALE < 10778 
EXSCALE will always be rounded off to one figure sig- 
nificance, If EXSCALE = 0, the X~scale will be computed 
by DRAW, This is called auto-scale, 
Y~scale in units per inch, with properties corresponing 
to those of EXSCALE, 
Distance, in inches, of the X-axis from the bottom of 
the graph (0 < IXUP S IHIGH)., This argument will be 
ignored unless MODEXAX = 2, see below. 
Distance, in inches, of the Y~axis from the left of the 
graph (O © IYRIGHT S IWIDE), This will be ignored unless 
MODEYAX = 2, see below. 
Determines the mode of the X-axis location: 
The X-axis will be located automatically by DRAW, with 


the origin of Y on the graph. 
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MODEYAX s 
IWIDE: 
IHIGH: 


IGRID: 


LAST: 


= 1 


=2 


The X-axis will be automatically located by DRAW, 

with the origin of Y removed (in one's imagimtion) an 
integer number of inches above or below the graph, if 
this is appropriate, This option can be used only if 
the Y~scaling is automatic (YSCALE = 0). ¥ 

The X~axis location will be as specified by IXUP, 
Determines the mode of Y-axis.location in the same way 
as MODEXAX, above, governs the X-axis location. 

Width of graph in inches (1< IWIDES 9). If IWIDE 

is out of this range, a value of 8 will be assumed, 
Height of graph in inches (1 S IHIGH £15). If IHIGH 
is out of this range, a value of 8 will be assumed. 

If IGRID =1, a 1" x 1" grid will be superimposed on 
the graph, This is useful only if plain paper is used 
on the CalComp Plotter, 


Indicates to the calling program whether the previous 


“plot was completed successfully. The codes are: 


Last plot was completed successfully. 
Last plot was not completed successfully. 
Last plot was not completed successfully, and no further 


graph output will be attempted until DRAW is next 


‘entered with MODCURV = 1 or 0, 


An attempt was made to enter DRAW with MODCURV #1 or 
O while the error lockout was set, 
This argument must always be a variable name and never 


&@ number in the call statement, 


Notes and Comments: 


a, The graph scales-and, if MODEXAX = 1 and/or MODEYAX = 1, the 
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amounts of origin offset - are always output as part of the graph title. 

b. Hach time a graph is completed, a message to this effect is 
printed on both the standard output and the console typewriter, 

c. There are internal checks of the input to DRAW to prevent incorrect 
use, If an input error is detected, an attempt will be made, where 
possible, to complete the plot. If an argument is "corrected" in this 
process, the user will be so informed on the standard output. If it is 
not possible to complete the plot, the user will be informed of the reason 
by a message on the standard output. 

d, If part or all of a curve would fall more than 0.6" laterally 
beyond the ends of the X-axis, or 0.7" vertically beyond the ends of -the 


Y-axis, the X and/or Y ordinates will be limited so that the curve will 


,, typically become a line along part or all of the boundary of the graph 


as here defined, 

e. If one or more points of a point plot would fall outside the graph 
area, the plot of that point, or poirfts will be inhibited... The number of 
such points will be reported to the user on the standard output. 

f. It should be pointed out that the X and Y scaling and axis locating 
processes are entirely independent, so that, for example, X might be 
auto-scaled, while the Y-scale is specified, At the same time the X-axis 
might be located automatically, while the Y-axis location is specified, 

ge It mist be remembered that the scales and axis locations of a 
multi-plot graph are set when DRAW is called for the first time (with 
MODCURV = 1). Thus the user mst attempt, at that time, to achieve 
scaling and axis location which will be appropriate to all the plots he 
intends to make on the one graph, Particularly if the automatic features 
of DRAW are selected, foresight will be demanded of the user in this 


respect. 
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h. Auto-Scale Properties: 

The scale factor is chosen from amongst the values 1,2, or 5 units 
per inch, or some power of 10 times one of those values. A curve, or 
set of points which is plotted with auto-scale will normally lie entirely 
within the graph area as defined in 3.d., above, The only exception may 
occur if an axis is placed, by the user, along one edge of the graph 
(e.g., IXUP = 0, MODEXAX = 2). In such a case, points "outside™ the axis 
are not considered in the selection of a scale factor (e.g., negative Xj 
do not affect the choice of scale when IXUP = 0). If automatic axis 
location as well as auto-scale is selected, the plot, if it does not fill 
the graph area, will be placed as far as possible towards the bottom-left 
of the graph area consistant with the fact that the axes can be set in 
increments of 1" only. 
5. Space Required: 3960 cells (excluding the input arrays). 
6. Temporary Storage: None. 
7. &rror Print-Outs: There are a large number of possible error print- 
outs. These are all self-explanatory. 
8. Error Returns: All error returns are preceded by self-explanatory 
error printouts. An error indication is transferred back to the calling 
program through the argument LAST, 
9. Error Stops: None. 
10. Tape Mountings: Logical Tape #8 will receive the binary graph output. 
The standard monitor output will receive the messages to the user, 
11. Output Format: The format of the binary graph output records on 
magnetic tape is described in reference 1. The only difference is that 
in this program the interpolation option is i Seheeea (Set to zero in 
the graph output record), See reference 2, 


12, Selective Jump and Stop Settings: None, 
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13, ‘ining: Variable, depending upon the number of points and the 
options chosen, Typically less than one second per curve or point plot. 
14. Accuracy: The accuracy of results is equal to the resolution of 
the CalComp Plotter, that is, 0.01" in both the X and Y directions, 

15. Equipment Configuration: CDC 1604 with FORTRAN 60 compiler and 
Library. A CDC 160 or 160A with CalComp 165 Plotter is needed for the 
off-line plotting using the appropriate GRAFPLOT program, 

De REFERENCES: 


1. Weir, Maurice D., Spritzer, Milton and McIlhenny, D.W., "160-A 
Graph Plot Program,® Ident *BOO1, SWAP Library, 15 August 1962. 


2. Hogg, R. L. and Glover, D.C., "160 Grafplot Routine*, Writeup 
Available from Computer Facility, U. S. Naval Postgraduate 
School, 1 April 1963. 


3. U.S. Naval Postgraduate School, Thesis, "Control System 


Programming, Remote Computing and Data Display,*® by Robert lee 
Hogg and Dennis C. Glover, 1963. 


N.B,. References 1 and 2 are included in Reference 3 as 
Enclosures 2 and 1, respectively. 
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APPENDIX E 
PLOTS OF DOCK TRIAL DATA 


Scale of stress has been omitted for reasons of security 
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